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TITULO 

COMPOSICI6N COHERENTE DE SENALES POR CORRECCI6N FOCAL 
PROGRESIVA. 

5 SECTOR DE LA TECNICA 
Tecnologfa ffsica 

Sistema electronico para Ecografia clmica y Ensayos No Destructivos. 
Otras aplicaciones: radar, sonar y seismologia 

1 0 ESTADO DE LA TECNICA 

En diversos ambitos de aplicacion, un conjunto de N transductores o array, 
recibe senales de una o multiples fuentes ubicadas a cierta distancia y 
orientacion relativa al array. Esta situacion es comun en radioastronomfa, 
radar, sonar, sismologia y, particularmente, en el ambito de la tecnologfa 
15 ultrasonica basada en arrays (ecograffa clmica o industrial, en este caso 
aplicada a los Ensayos No Destructivos), donde N se encuentra, tipicamente, 
en el rango 32 a 256 o mas. 

En estas situaciones, cada uno de los N elementos del array recibe la senal 
generada por cada fuente con cierto retraso temporal respecto a los restantes 

20 elementos. La funcion de un conformador de haz es obtener una unica senal a 
partir de las N recibidas en la que los retrasos temporales para cada posible 
distancia en una orientacion determinada se compensen, aumentando la 
capacidad de discriminacion del sistema (mayor resolucion, relacion senal/ruido 
y contraste). Repitiendo este proceso en una diversidad de direcciones 

25 mediante barridos lineales o angulares, se obtiene informacion suficiente para 
poder componer una imagen 2D o 3D de las fuentes de senal. En el caso 
particular de los ultrasonidos, las N senales recibidas son ecos producidos por 
discontinuidades en el medio que ha sido insonificado por un emisor, que 
puede ser el mismo array receptor (o un subconjunto del mismo) u otro 

30 transductor (array o monoelemento). 
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En recepcion, es deseable compensar los retardos en todos los puntos o focos 
a lo largo de la Imea de propagacion en una direccion determinada. Esto 
requiere modificar continuamente los retardos en el tiempo, para adaptarlos a 
la llegada de los ecos desde distancias dinamicamente variables a cada 
5 elemento. Esta tecnica de focalizacion dinamica, proporciona la mayor 
resolucion, contraste y relacion senal/ruido. 

Siendo tecnicamente imposible modificar en forma continua los retardos, es 
esencial que su discretizacion se realice con una excelente resolucion temporal 
para cada elemento del array y para cada distancia, tfpicamente del orden de 
10 1/16 a 1/64 el periodo fundamental de la serial. 

Este hecho ha motivado la busqueda de metodos para la focalizacion dinamica 
en recepcion sencillos y eficientes que, a lo largo de los anos, se han ido 
adaptando a las disponibilidades tecnologicas. Asi, se han descrito tecnicas 
basadas en dispositivos CCD accionados por relojes cuyo periodo determina el 
15 retardo [Macowski, "Ultrasonic electronic lens with reduced delay range", US. 
Pat. Nr. 4.058.003, Nov. 1977] o en Imeas de retardo conmutables [Maslak, US. 
Pat. Nr. 4.550.607, "Phased array acoustic imaging system", Nov. 1985]. Estos 
metodos introducen artefactos de conmutacion, son voluminosos y 
proporcionan una resolucion limitada estando, practicamente, en desuso. 

20 Otra alternativa, asimismo analogica, se basa en la utilizacion de mezcladores 
y desplazadores de fase que realizan los retrasos temporales en el dominio de 
la frecuencia [Maslak et al., "Dynamic focused linear phased array acoustic 
imaging system", US. Pat. 4.699.009, Oct. 1987]. Uno de los problemas 
asociados a este tipo de tecnicas es la dificultad de su adaptacion a 

25 transductores de distintas frecuencias. 

Por otra parte, se ha propuesto la utilizacion de rotadores de fase, basados en 
el algoritmo Cordic [R. Andraka, "A survey of CORDIC algorithms for FPGAs", 
Proc. 1998 ACM/SIGDA VI Int'l Symposium on FPGAs, Monterrey, CA, pp. 191- 
200, Feb. 1998]), y su aplicacion a la focalizacion dinamica se encuentra en [O' 
30 Donnell et al., "Method and apparatus for digital phased array imaging", US Pat. 
Nr. 4.983.970, Ene. 1991] y en [Engeler et al. "Dynamic phase focus for 
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coherent imaging beam formation", US Pat. 5.111.695, May 1992]. El hardware 
involucrado es significative 

La disponibilidad de tecnologia VLSI permitio nuevas propuestas basadas en 
retardar las senales una vez digitalizadas. Cada elemento del array se asocia a 
5 un canal de adquisicion y procesamiento (focalizacion) independiente. Las 
salidas de los N canales se suman para obtener la senal resultante (suma 
coherente). Para operar en tiempo real, cada canal dispone de dos espacios de 
memoria: uno para datos, donde se almacenan las muestras adquiridas, y otro 
de focalizacion que establece los retardos a aplicar. 

10 En un principio, la resolucion de los retardos se limito al propio periodo de 
muestreo de las senales, determinado por el criterio de Nyquist, alrededor de 
1/3 a 1/5 del periodo de la senal, utilizando registros de desplazamiento de 
longitud variable [Welles II et al., "Architecture for ultrasonic imaging", US Pat. 
4.796.236, Ene 1989]. La baja resolucion temporal produce lobulos de 

15 cuantificacion de retardos que reducen el rango dinamico, contraste y relacion 
senal/ruido de la imagen [D. K. Peterson, G. S. Kino, "Real-time Digital Image 
Resconstruction: A Description of Imaging Hardware and an Analysis of 
Quantization Errors", IEEE Trans. Sonics Ultrasonics, 31, pp. 337-351, Jul. 
1984]. 

20 La opcion de utilizar una mayor frecuencia de muestreo no es aconsejable por 
el mayor coste de los componentes, consumo energetico y volumen de datos a 
almacenar y procesar. Para solventar este problema se han propuesto diversas 
tecnicas que, basicamente, pueden agruparse en dos: 

1) Realizacion de un retardo grueso con una resolucion temporal equivalente al 
25 periodo de muestreo, y otro fino mediante interpolacion de la senal recibida 
entre dos muestras consecutivas, con una resolucion, tipicamente, entre 1/4 y 
1/16 del periodo de muestreo. 

A este grupo pertenecen las patentes: [Odell, "Digital beamforming and filtering 
circuit", U.S. Pat. Nr. 5,268,877, Die. 1993], [D.Lipschutz, "Time multiplexed 
30 digital ultrasound beamformer"US.Pat. 5,469,851, Nov. 1995], [Beaudin et al. 
"Delay generator for phased array ultrasound beamformer", US. Pat. 5.522.391, 
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Jun. 96], [Wright et al. "Method and apparatus for receive beamforming 
system", US Pat. 6.042.547, Mar. 2000] y [M.H. Bae, "Ultrasonic signal focusing 
method and apparatus for ultrasonic imaging system", US Pat. 6.231 .51 1 , May. 
2001]. 

5 Esta tecnica, con diversas variantes para la realizacion fisica de los retardos 
grueso y fino, incluida la de rotacion de fase mediante algoritmos CORDIC, 
requiere invariablemente una etapa de interpolacion y un selector de muestras 
interpoladas por cada canal, con la consiguiente complejidad y errores 
asociados. 

10 2) Adquisicion directa de las senales en los instantes que corresponden a la 
llegada de la senal de cada foco a cada elemento. Para ello, se controla la fase 
del reloj de muestreo de cada canal con una resolucion temporal equivalente 
entre 1/16 y 1/64 el periodo fundamental de la senal. Este metodo hace 
innecesaria la etapa de interpolacion, eliminando la electronica y las fuentes de 

15 error correspondientes. Las muestras se pueden adquirir sobre una memoria 
FIFO, de la que se extraen para realizar la suma coherente [Chang et 
al., "Bandwidth sampling technique for digital focusing in array imaging 
systems", US. PAT. 5,581,036, Die. 1996]. 

El mayor problema asociado a esta tecnica es la realizacion de los generadores 
20 de reloj de muestreo de fase controlada con la suficiente precision. En 
[O'Donnell et al., "Method and apparatus for fully digital beam formation in a 
phased array coherent imaging system", US Pat. Nr. 4.809.184, Feb. 1989] y 
[Corl, "Digital ultrasound system with dynamic focus", U.S. Pat. 4,974,211, Nov. 
1990], se generan m relojes derivados de un reloj maestro de frecuencia m 
25 veces superior a la de muestreo, cada uno de ellos con un desfase 2n/m; un 
selector del reloj por cada canal, accionado por el contenido de la memoria de 
focalizacion, permite asociar en cada instante uno de los relojes al de 
muestreo. La distribucion de multiples relojes y su seleccion con coherencia 
temporal es un problema, especialmente para las mayores resoluciones 
30 temporales. 

Otra posibilidad es almacenar el reloj de muestreo de cada canal en una tabla 
como una secuencia de 1s y 0s, como se ha sugerido en [J.H. Kim, T. Song, 
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S.B. Park, "Pipelined Sampled-Delay focusing in ultrasound imaging 
systems, Ultrasonic Imaging 9, pp.75-91, 1987], pero no es practica por la 
cantidad de memoria necesaria para almacenar todos los relojes de muestreo 
de todos los canales para multiples Imeas de barrido. 

5 Siguiendo otra alternativa, se calculan en tiempo real los instantes de muestreo 
mediante circuitos especializados asociados a cada canal. Esta tecnica se 
describe en [Park et al., "Real time digital reception focusing method and 
apparatus adopting the same", US Pat. 5.669.384, Sep. 1997], [Bae, " Focusing 
delay calculation method for real-time digital focusing and apparatus adopting 
10 the same" US Pat. 5,386.881, Feb. 1998] y [Petrofsky, "Method and apparatus 
for distributed focus control with slope tracking", US Pat. 5.724.972, Mar. 1998]. 
Todas, con diferentes variantes, se basan en el algoritmo del punto medio 
descrito en [Bae et al., An Efficient Real Time Focusing Delay Calculation in 
Ultrasonic Imaging Systems, Ultrasonic Imaging, 16, pp. 231-248, 1994]. 

15 Aunque se trata de una tecnica interesante por sus posibilidades de integracion 
en un ASIC, presenta diversos inconvenientes entre los que destacan: la 
necesidad de programar un gran numero de parametros para cada Imea de 
barrido y cada canal, la complejidad del hardware para obtener una buena 
resolucion, menor eficiencia de utilizacion de area de silicio frente al requerido 

20 por sistemas basados en memoria por la baja regularidad del hardware y, 
finalmente, requiere un compromiso entre obtener una alta velocidad de calculo 
o una reduccion del rango de aplicacion. 

El metodo y sistema electronico objeto de la presente invencion pertenece a 
este grupo y resuelve los problemas descritos mediante una tecnica mixta, que 
25 combina una memoria de focalizacion de reducido tamano, en la que solo se 
requiere una fraccion de bit por cada foco y canal, con el calculo en tiempo real 
del instante de muestreo con elevada resolucion mediante circuitos sencillos y 
regulares de facil integracion VLSI. 



30 
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DESCRIPCI6N DE LA INVENCI6N 
INTRODUCCl6lM: 

Se describe un conformador de haz que realiza la composicion coherente de 
las senates procedentes de una direccion de propagacion determinada, 
5 recibidas por un conjunto de N transductores receptores o array. La 
metodologia y aparatos descritos son aplicables en campos como imagen 
ultrasonica (ecografia y Ensayos No Destructives), radar, sonar, acustica, 
seismologfa y otros. 

El conformador de haz aqui descrito opera con arrays cuyos elementos pueden 
10 estar distribuidos arbitrariamente. En particular, realiza la focalizacion dinamica 
con deflexion del haz en acimut y elevacion utilizando arrays bidimensionales. 

Basa su operacion en muestrear la serial recibida por cada elemento del array 
en los instantes que corresponden a los de llegada de la senal desde multiples 
focos situados a lo largo de una direccion de propagacion dada. Dentro de los 

15 Ifmites establecidos por los principios que guian la presente invencion, el 
instante de adquisicion de la senal procedente de cada foco y recibida por cada 
elemento del array se determina con un error absoluto, respecto al valor 
exacto, inferior a medio periodo de un reloj maestro. La frecuencia del reloj 
maestro puede elegirse para reducir la incidencia negativa de los lobulos de 

20 cuantificacion de los retardos en la relacion senal/ruido, contraste y rango 
dinamico de la senal resultante. 

La frecuencia de muestreo es independiente de la de la senal, siempre que se 
verifique el criterio de Nyquist, tanto en radiofrecuencia como en banda base, y 
de la frecuencia del reloj maestro. Con tecnicas de sobremuestreo (conversion 
25 AL), la frecuencia de digitalizacion puede ser la mas baja compatible con la 
relacion senal/ruido deseada. En particular, puede utilizarse el reloj maestro 
modificado segun el metodo descrito para realizar directamente el muestreo 
AL. 

Al reducir la frecuencia de muestreo a la mas baja posible, compatible con el 
30 ancho de banda de las senales y/o relacion senal/ruido, se reducen el coste, 
consumo energetico y volumen de datos a procesar. En particular, con el 
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metodo que constituye una parte esencial de la presente invencion, solo se 
adquieren las muestras necesarias para realizar la suma coherente, sin que 
sea necesario un proceso de interpolation o la adquisicion de informacion 
redundante. 

5 Ademas, el metodo de conformacion que se describe y que constituye otra 
parte esencial de la presente invencion, posee unos bajos requisites de 
memoria de focalizacion, que pueden llegar a ser de una fraccion de bit por 
cada muestra adquirida en cada canal. Esta alta eficiencia facilita que el 
metodo pueda aplicarse convenientemente a una diversidad de modos de 
10 operacion, como son el barrido lineal o sectorial, con aperturas densas o 
dispersas, y otras aplicaciones en las que se requieren multiples descriptores 
de la focalizacion dinamica para cada configuracion de la apertura activa del 
array. 

Dichas caracteristicas posibilitan la integracion conjunta de la memoria de 
15 focalizacion, la de adquisicion de datos, la logica de generacion de los relojes 
de muestreo y la de otras funciones auxiliares como apodizacion, apertura 
dinamica y control, de multiples canales, en un unico dispositivo VLSI o 
estandar configurable tipo FPGA, con la consiguiente reduccion de coste sin 
perdidas de calidad. 

20 

Descripcion del metodo. 

Por sencillez, la exposicion se realiza para un array unidimensional de N 
elementos, no necesariamente equiespaciados, considerando una aplicacion 
de ultrasonidos utilizando la tecnica de pulso-eco en un medio homogeneo con 
25 velocidad de propagacion c. El metodo que se describe a continuacion se 
denominara genericamente tecnica de Correction Focal Progresiva o CFP, que 
describe la forma en que opera. 

La Figura 1 muestra un diagrama de la geometria general con los 
condicionantes anteriores. El array, de N elementos, numerados 1, 2, 3, k, 
30 AM y N, se situa sobre el eje de las x, con el origen de coordenadas 

arbitrariamente definido, por comodidad, en el centra del array. Se desea hacer 
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una adquisicion de senal en la direccion angular 0 referente al eje z. Para ello, 
se definen multiples focos, a intervalos AR en la direccion de propagacion 
(llegada de la serial). El metodo propuesto combina las senales procedentes de 
cada foco y recibidas por los diferentes elementos del array de forma que se 
5 obtiene una salida que representa la senal que seria proporcionada por un 
transductor del tamano del array focalizado en todos y cada uno de los focos. A 
este proceso se le denomina composition coherente de las senales recibidas y, 
al resultado, senal focalizada dinamicamente con deflexion. 

Los intervalos AR se eligen convenientemente en funcion del periodo de un 
10 reloj maestro T x =1/fxy de la velocidad de propagacion c como: 

AR=vcT x /2 (1) 

donde v es un numero entero que expresa el numero de periodos del reloj 
maestro que hay entre dos focos consecutivos en pulso-eco. Por otra parte, los 
intervalos AR deben ser lo suficientemente pequenos como para que se 

15 verifique el criterio de Nyquist: si T s es el intervalo temporal entre muestras, y 
T R es el periodo fundamental de la senal, el criterio de Nyquist se verifica, 
tfpicamente, para T s < Tr/3, considerando el ancho de banda de las senales 
habituales y la adquisicion en radiofrecuencia. Preferentemente, el intervalo 
temporal entre muestras es un multiplo del periodo del reloj maestro, esto es, 

20 T s = sT x , con s entero. 

En general, puede adquirirse un numero m de muestras entre dos focos 
consecutivos, eligiendose preferentemente el valor v en funcion del intervalo 
entre muestras Ts = sT x para que v = ms, con my s enteros. En una 
realizacion particular en la que m=1, cada muestra coincide con un foco, de 

25 modo que todas las muestras se focalizan dinamicamente. Cuando m>1, 
estrictamente solo se realiza la focalizacion dinamica en los focos, aunque la 
profundidad de foco permite, en general, suponer bien focalizadas las muestras 
intermedias. Una variante del metodo, denominada tecnica de Correction Focal 
Progresiva Variable o CFPV, varfa dinamicamente el numero de muestras 

30 adquiridas entre dos focos, reforzando la idea de mantener una buena 
focalizacion dinamica en todas las muestras adquiridas. 
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Por otra parte y, como es bien conocido, el valor jll = Tr/T x o relacion entre el 
periodo de la serial recibida y la resolucion temporal del sistema, afecta al 
rango dinamico y contraste de las imageries construidas con multiples haces 
focalizados. En general, se considera aceptable un valor \y >16 (ver, por 
5 ejemplo, [D. K. Peterson, G. S. Kino, "Real-time Digital Image Resconstruction: 
A Description of Imaging Hardware and an Analysis of Quantization Errors", 
IEEE Trans. Sonics Ultrasonics, 31, pp. 337-351, Jul. 1984]). 

Es importante que la resolucion temporal determinada por T x sea 
independiente del periodo de muestreo T s que se elige en base al criterio de 

10 Nyquist y, por tanto, se relaciona con el periodo de las senales recibidas, T R . 
Cuando el muestreo se realiza en radiofrecuencia se elige, en general, Ts » 
Tr/3o T s ~ Tr/4 para reducir el flujo de datos a procesar y el consumo 
energetico, esto es, la frecuencia de muestreo mas baja compatible con el 
criterio de Nyquist. En tal caso, con s> 6 o s > 4, respectivamente, se tiene 

15 que (lx > 16. Una particularidad de las tecnicas CFP y CFPV es que permiten 
operar a la mayor parte de la electronica a la frecuencia de muestreo, 
reservandose la mas elevada del reloj maestro para muy pocos elementos que 
realizan el ajuste del instante de muestreo. 

Para un foco F, situado en la coordenada polar generica (R, Q) y, para el 
20 elemento k del array situado en la coordenada cartesiana (x k , 0), por el teorema 
del coseno (Figura 1 ): 



donde L k (R,Q) representa la longitud del camino recorrido por la serial desde el 
foco F situado en (R,Q) al elemento k. El tiempo de vuelo de un pulso 
25 ultrasonico desde su emision arbitrariamente elegido en el origen de 
coordenadas al foco F y de este al elemento k, es: 



En el caso en que el origen de tiempos se situe en una coordenada diferente 
(por ejemplo, cuando el espacio inspeccionado sea iluminado desde un emisor 




(2) 




(3) 
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diferente del propio transductor receptor), el primer termino de la Ecuacion (3) 
debera ser modificado para representar el tiempo de vuelo desde el emisor al 
foco F. 

En general, se denomina 7V/ al tiempo invertido por la senal correspondiente al 
5 foco Ft en alcanzar el elemento k del array. Esta notacion generica se refiere, 
tanto a la geometria indicada en la Figura 1, como para otras situaciones en las 
que los elementos del array se distribuyen de forma arbitraria, el medio de 
propagacion no es homogeneo o los focos se situan en otras posiciones. No 
obstante, para facilitar la exposicion, se mantiene la geometria de la Figura 1, 
10 el medio homogeneo con velocidad de propagacion c y las posiciones de los 
focos se alinean formando un angulo G respecto al array situandose, de 
momento, a intervalos regulares AR a partir de una distancia inicial R 0 en la que 
se coloca el primer foco F 0 . 

La variacion en el instante de llegada al elemento k de sehales procedentes de 
15 dos focos consecutivos F h1 y F, situados a distancias R r AR y R f del origen de 
coordenadas, respectivamente, esta dada por la diferencia: 

AT U (R. ,e) = T k (R r e)-T k (R . -AR,Q) (4) 
que, con AR suficientemente pequeno, es: 



AT k .(R.,Q) = AR 



dT k (R,Q) 



dR 



(5) 



R = R. 
i 



20 Derivando la Ecuacion (3) respecto a R y sustituyendo el valor de AR dado por 
la Ecuacion (1): 



1 + 



R. - x 1 senQ 
i k 



vT 



Ft 



2 +x\ -IR.x.senQ 



X 



(6) 



Esta es una funcion monotona creciente entre R=0 y R— » «> e independiente de 
la velocidad de propagacion, con los siguientes limites: 
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f 




VT 



X 



lim I 



1- 



senQ 



2 



>0 



(7) 



v 



J 




X 



(8) 



donde, por claridad de notacion, AT k = AT k (R,Q). El valor x^\x k \ de la Ecuacion 
(7) representa el signo de x k , mientras que 0 interviene con su signo. En 
5 cualquier caso, el menor valor posible para AT k es cero, independientemente 
de la posicion del elemento k. 

De este modo, los intervales entre los instantes de muestreo que corresponden 
a dos focos consecutivos estan limitados por los valores expresados en las 
Ecuaciones (7) y (8). En particular, para el caso 9=0 (no se realiza deflexion), la 
1 0 funcion crece monotonamente entre vT x /2 yvT x cuando la distancia a la que se 
situan los focos varfa de 0 a oo. Cuando se reciben las senales con deflexion 
del haz, el rango de variacion es mayor y, en un caso extremo de deflexion a 
90°, entre 0yv7 x . 

En una posible realizacion, las posiciones de los focos quedan definidas por 
15 valores comprendidos en el intervalo [0, v], que representan los valores 
incrementales dados por la Ecuacion (6) en periodos de reloj maestro. En este 
sentido, si el numero entero v puede expresarse con b bits, cada retardo 
incremental para cada foco y para cada elemento puede representarse con 
esta misma anchura, independientemente de la posicion ocupada por el 
20 elemento del array, de la distancia a la que se situe el foco y del angulo de 
deflexion. El error absolute con el que se determina el instante de muestreo es 
inferior a medio periodo de reloj maestro al redondear al entero mas proximo el 
resultado de AT ki (R,Q)/T x dado por la Ecuacion (6). El inconveniente de esta 
opcion es que requiere un elevado numero de bits para expresar los retardos 
25 incrementales para cada foco y elemento, con las consiguientes demandas de 
capacidad de memoria. 
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Inversamente, si se dispone de b bits para representar cada retardo 
incremental, el rango de valores posibles puede elegirse entre dos situaciones 
extremas: 



(v-2" + iy x <AT k <vT x 
0< AT k <{2 b ~W X 



(9a) 



(9b) 



La forma (9a) resulta util para realizar la focalizacion dinamica a partir de una 
distancia minima Rmin, esto es, para R > Rmin, mientras que la forma (9b) seria 
preferible para focalizar dinamicamente en las regiones mas proximas al array, 
esto es, para R < R m ax- Entre uno y otro caso pueden darse rangos intermedios 
mas adaptados a una aplicacion particular, donde R min <R< R ma x- 

Sin embargo, con la Tecnica de Correccion Focal Progresiva se prefiere la 
forma (9a) que no limita la maxima distancia de validez. Para las regiones mas 
proximas al array en las que no es posible verificar la desigualdad de la 
izquierda de la Ecuacion (9a), la tecnica propuesta proporciona una serie de 
alternativas que seran descritas con detalle mas adelante. Asi, pues, llamando 
en la Ecuacion (9a): 



a = 2 b -1 



puede escribirse en forma mas compacta: 

(v - a)T x < AT k <vT x 



(10) 



(11) 



La Ecuacion (6) siempre verifica la desigualdad de la derecha de la Ecuacion 
(11), de acuerdo con la Ecuacion (8). Para que tambien verifique la de la 
izquierda: 



1 + 



R - jc t send 
k 



ft 



2 + x| - 2Rx k senQ 



vT 

JL>(v-a)T x 



(12) 



Con: 
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0 9 

xf cos z e 



2 + x} - 2Rx, senQ = (r - x^send )Jl + — - ~ 

[R-x,senQ) 2 



la Ecuacion (12) queda como: 



1 + 



1 



i + 



2 2 
x, cos 0 

k 

(r - x k senQ j 2 



V ^ 

— >v - a 
2 



donde, haciendo los cambios de variable: 



r = R — x , senQ 
k 



sb obtiene: 



(13) 



(14) 



(15) 



10 



15 



1 + 



1 



— >v - a 
2 



Despejando rde la ecuacion anterior se obtiene: 

z k (y -2a) 



r> 



v >2a 



(16) 



2-Jav - a 2 

que, deshaciendo los cambios de variable de las Ecuaciones (15) conduce a: 



(17) 



R>R Q (x k ,e) = 



v -2a 



l4av - a 2 



x 7 cos9 
k 



x, send 
k 



v >2a 



(18) 



donde |-| representa el valor absoluto y considera el caso mas desfavorable de 
signos relativos de x k y 9. La ecuacion (18) representa una cota superior para 
la distancia minima Ro(x k , Q) a partir de la cual los intervalos temporales para 
muestrear la serial recibida por el elemento k con un angulo de deflexion 9 
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pueden representarse con un numero entero de b bits, con un error inferior a 
medio periodo T x de un reloj maestro. 

Si la Ecuacion (18) se verifica para todos los elementos del array a todas las 
distancias y angulos de la region a inspeccionar, el retardo inicial para adquirir 
5 la primera muestra, correspondiente al foco F 0 , en el canal /cviene dado por: 

T oW>=i{ R o + K +x l - 2R o x k senQ ) (19) 

donde R 0 = R ofa, Q) dado por la Ecuacion (18). Obviamente la primera 
adquisicion puede realizarse a una distancia Ri > R 0 , en cuyo caso el valor de 
To(x k , Q) se calculara para este nuevo valor. 

10 A partir de este instante, los intervalos de muestreo entre focos presentan un 
error acotado a medio periodo del reloj maestro, pudiendo codificarse para 
cada foco F/de coordenadas (R h Q) en el canal /ccomo: 



T 



(20) 



donde [ ]u representa la funcion de redondeo al entero mas proximo y AT k} (R h 
15 Q) esta dado por el valor de la Ecuacion (6). Sustituyendo en (20) las 
Ecuaciones (1 0) y (1 1 ): 

0<e^<2 fe -l (21) 

que indica que los intervalos de muestreo pueden codificarse con b bits, 
produciendo un error inferior a medio periodo de reloj maestro. El valor Q ki es el 
20 codigo de correction focal que determina el numero de periodos de reloj 
maestro que debe adelantarse la adquisicion de la serial correspondiente a I 
foco / en el canal k. 

La adquisicion de un numero m > 1 de muestras entre focos aumenta la 
eficiencia de memoria al compartir un unico codigo de correccion focal de b bits 
25 entre varias muestras, con una densidad de b/m bits por muestra. Si, por 
ejemplo, 6=1 y m=4, se requieren 0'25 bits para focalizar cada muestra 
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adquirida. Una buena eleccion para reducir el error en los intervalos de 
muestreo consiste en distribuir uniformemente el adelanto expresado por el 
codigo de correccion focal Q k i entre las m mue stras. Para esto, en una 
realizacion preferida, se calculan los valores: 



v-1 
m 



\<j<m (22) 

n 



donde q=Qik representa, genericamente, el codigo de correccion focal para 
simplificar la notacion y []u es la funcion de redondeo. El valor de u(j, q) 
representa el entero mas proximo del numero de periodos del reloj maestro 
entre el primer foco y la muestra ycon j=1, 2, m, consiguiendo una 
10 distribucion lo mas uniforme posible del adelanto q entre las m muestras. 

El intervalo temporal entre un foco cualquiera con la muestra que le sigue en el 
orden j=m es u(m, q) = v-q, esto es, la muestra j=m se adquiere exactamente 
en la posicion del siguiente foco, con un adelanto de q ciclos del reloj maestro 
respecto del intervalo nominal v entre focos, de acuerdo con lo expresado por 
15 el codigo de correccion focal. 

El intervalo entre muestras, en periodos del reloj maestro, viene dado por: 

Aw(l, q) = q) j — 1, por definicion 

Au(j, q) = u(j, q) ~ u(J -\,q) 2<j<m < 23 > 

En el caso mas habitual en el que el numero de muestras entre focos es mayor 
o igual al maximo adelanto a aplicar entre focos, esto es, q < 2*^-1 < m, al 

20 hacer una distribucion uniforme de q ciclos entre m muestras, el adelanto g, 
asignado a la muestra j solo puede tomar los valores 0 6 1 respecto al intervalo 
nominal v/m que se da para q=0. Este hecho facilita la realizacion de la 
distribucion uniforme del adelanto q entre un numero de muestras m > q 
mediante una tabla en la que se codifica con 1 bit el valor g(j,q), previamente 

25 calculado segun: 

gU,q)=v /m-AuU,q) l<j<m 9 0<q<2 b -1 (24) 
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Para clarificar el metodo, considerese un ejemplo en el que m=3, v=12=4 m, 
b=2. Los posibles valores de q son 0, 1 , 2 y 3 por estar codificado con 2 bits, y 
v/m=4. La Tabla I muestra los resultados que se obtienen al aplicar, con estos 
valores, las Ecuaciones (22), (23) y (24). Se observa que el adelanto gse 
5 distribuye uniformemente entre las muestras dentro de la resolucion de 1 
periodo de reloj maestro para todos los codigos de correccion focal. La 
densidad de la codificacion en este ejemplo es de 2/3 = 0.67 bits/muestra. 

Tabla I. 

10 Ejemplo de distribucion uniforme de q adelantos entre m muestras 





v-q 


u(1,q) 


u(2,q) 


u(3,q) 


Au(1,q) 


Au(2,q) 


Au(3,q) 


9(1 -q) 


9(2,q) 


9(3,q) 


q=0 


12 


4 


8 


12 


4 


4 


4 


0 


0 


0 


q=1 


11 


4 


7 


11 


4 


3 


4 


0 


1 


0 


q=2 


10 


3 


7 


10 


3 


4 


3 


1 


0 


1 


q=3 


9 


3 


6 


9 


3 


3 


3 


1 


1 


1 



Para cada valor de ghay m valores posibles g(j,q)> cada uno de 1 bit. Por 
consiguiente, todos los adelantos individuales pueden codificarse en una 
memoria de m-^ palabras de 1 bit (en el ejemplo anterior, basta una memoria 
15 de 12x1 bit). 

Para valores R< Ro(x k , Q) el error absoluto del intervalo entre instantes de 
muestreo puede superar el valor T x /2. En este caso, las tecnicas CFP y CFPV 
pueden operar de varias formas, destacando las siguientes estrategias: 

1. Activar una funcion de apertura dinamica en la que el elemento kno 
20 interviene hasta que se reciban senales desde una distancia igual o superior 
a Rofok, Q)- Una de las posibles realizaciones, produce para este elemento 
un numero de muestras nulas M z (k) equivalente a las que se producen a la 
salida del conformador entre el comienzo de la adquisicion T ini y el instante 
T 0 (x kf Q) dado por la Ecuacion (19): 
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T 0 (x k ,B)-T. 



ini 



T 



(25) 

T 



donde [ ]r representa la funcion de redondeo por exceso a un numero 
entero. 

Debe destacarse que la relacion D/max(R 0 ), donde D es el tamano de la 
5 apertura, determina el mini mo valor de F# que puede aplicarse con las 

tecnicas propuestas. En un ejemplo practico se muestra que F# min =1, 0'6 y 
0'4 para valores de b= 1, 2 y 3, respectivamente, muy competitivos con los 
proporcionados por otras tecnicas. 

2. Anticipar la adquisicion de senal en el elemento k a una distancia Ri < 
10 Rofok, Q) elegida de forma tal que se adquiera una muestra en R 0 (x k , Q)- Los 

instantes de muestreo anteriores a este tendran un cierto error superior al 
acotado, pero su efecto en la imagen puede ser limitado. Para ello, el 
instante de adquisicion de la primera muestra dado por la Ecuacion (19) 
queda modificado a: 

15 r Q M {x k ,e ) = t q (x k ? e ) - m z (k)T sl (26) 

donde se expresa que, para el elemento /c, se adquieren M z (k) muestras 
anteriores a la recepcion con error acotado y T S i es el periodo de muestreo 
utilizado, constante en esta region. El valor de T S i puede ser igual al 
periodo de muestreo de salida del conformador T s o, preferiblemente, un 
20 valor inferior para reducir en lo posible los errores en el instante de 

muestreo. 

3. Operar de forma mixta, donde un primer rango de distancias [R 2 , Ri) se 
cubre con apertura dinamica, otro [Ri, R 0 ) con un error limitado como en el 
punto anterior y, finalmente, un tercero [R 0 , ~>) en el que los errores de 

25 muestreo estan acotados dentro de medio periodo de reloj maestro 
aplicando la correccion focal progresiva. 

En una realizacion preferida que, ademas resulta de aplicacion general cuando 
la geometria del array es diferente, el medio no es homogeneo (por ejemplo, 
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cambios de velocidades de propagacion) o los focos se distribuyen a intervalos 
regulares pero siguiendo una trayectoria diferente a la considerada, se utiliza 
un algoritmo de calculo basado en los siguientes pasos: 

ALGORITMO 1 : 

5 PASOI: Se calculan los tiempos de propagacion T ki de la senal 
correspondiente a cada foco F h #= 0, 7, 2, hasta cada elemento k 
del array. El valor de T ki se obtiene en cada caso al considerar la 
geometria del sistema y la(s) velocidad(es) de propagacion. 

PASO 2: Para cada elemento k, se calculan los intervalos U ki entre instantes 
10 de llegada de la senal correspondiente a dos focos consecutivos, 

expresados en periodos del reloj maestro, como: 



T -T 
ki k, i - 1 

T 

. x j n 



i>\ (27) 



PASO 3: El valor v o intervalo nominal en periodos de reloj maestro entre 
focos se elige de los resultados de la ecuacion anterior como: 

15 v>max(C/ fe .) Vk,i (28) 

PASO 4: Si el numero de bits disponibles para codificar las correcciones 
focales es b, determinar para cada canal k el indice i=h k a partir del 
cual: 

v-U M <2 b -\ i>h k (29) 

20 PASOS: Para los focos F/Con /< h k , el numero de muestras M z (k), para 
operar con apertura dinamica o adquisicion anticipada en el canal k, 
es: 

M z {kY{h k -X)m (30) 

PASO 6: Los codigos de correccion focal para cada elemento k 
25 correspondientes al foco / > h k se calculan como: 

a,=v-^« *>K ( 31 ) 
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10 



15 



que expresan el numero de ciclos de reloj maestro que debe 
adelantarse la adquisicion de la serial correspondiente al foco / en el 
canal k. 

En otra realization preferida, durante la adquisicion de serial, la distancia entre 
dos focos consecutivos F h1 y F, se varfa en momentos seleccionados 
incrementando el numero de muestras que se adquieren entre ellos. Esta 
variante se denomina Teenies de Correction Focal Progresiva Variable o 
CFPV, en la que el numero de muestras adquiridas entre dos focos 
consecutivos pasa de ser una constante m a una funcion del numero de foco 
actual, esto es M=m(i). 

La Ecuacion 18 indica que el valor de R 0 crece con el de v. Por su parte, el 
valor minimo de v esta determinado por la Ecuacion 1 cuando se adquiere una 
muestra por foco, con un intervalo temporal entre muestras T s = sT x . Iniciando, 
pues la adquisicion con M=m(1)=1 a partir de una distancia Ro(V determinada 
por el valor v = s a la que, para todos los elementos: 



Rod)* 



s -2a 



l4as-a 2 



x 7 COS0 
k 



+ 



x 7 send 
k 



VJfc 



(32) 



donde a = 2 b - 7 . La adquisicion continua manteniendo M=m(1)=1 hasta 
alcanzar cierto foco a una distancia R=R 0 (2) a partir de la cual se verifica, para 
todos los elementos: 



20 



R 0 (2)> 



2s -2a 



l^2as-a 2 



x 7 COS0 
k 



+ 



x 7 send 
k 



Vk 



(33) 



25 



donde comienza un tramo en el que pueden adquirirse 2 muestras por foco, 
esto es, M=m(2)=2, manteniendo el error de muestreo acotado. Cierto tiempo 
despues, se alcanzara otro foco en R=R 0 (3) a partir del cual puede operarse 
con M=m(3)=3 y asf sucesivamente, de forma que, genericamente, a partir de 
un foco situado en R=R 0 (m) a partir del cual pueden adquirirse m muestras, 
verificandose: 
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ms - 2a 



x, COS0 



+ x 7 send 



(34) 



2i 



ams - 




En definitiva, con la Tecnica de Correccion Focal Progresiva Variable el numero 
de muestras que se adquieren entre focos se incrementa con el tiempo, de 
forma que los codigos de correccion focal puedan seguir expresandose con b 
5 bits. Debe destacarse que la condicion (34) debe ser verificada por todos los 
elementos para cada valor de m. 

En una realizacion preferida, se asocia a cada foco F, un campo de 1 bit, 
denominado codigo de alejamiento focal J h Cuando JpO se mantiene el actual 
valor m(i)=m(i-1) y, cuando J/=7 se incrementa: m(i)=m(i-1)+1. El algoritmo 2 
10 proporciona un posible metodo para el calculo de los valores de J,: 

ALGORITMO 2: 

PASO 1: En funcion de la distancia minima establecida por la aplicacion, R min , 
determinar el mayor valor m=m 0 para el cual se verifica la Ecuacion 



PASO 2: Inicializacion: m=m 0 . Variables auxiliares: /=7 (numero actual de 
foco), i=1 (numero actual de muestra). 

PASO 3: Repetir: 

a) Calcular Ro(m) dado por la Ecuacion (34) y el valor AR = cTs/2, 



15 



(34). 



20 



siendo c= velocidad de propagacion, T s = periodo de 



muestreo. 



b) Si Rmm + (i+m) AR > (m+1)R 0 



{ J f = 1;m=m+1;} 



en otro caso { J f = 0;} 



25 



c) f=f+1; 



d) i=i+m; 



hasta que / > N s , siendo N s = numero de muestras a adquirir. 
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El numero de codigos de alejamiento focal iguala al de focos para cualquier 
adquisicion, por lo que puede realizarse facilmente asociando a cada codigo de 
correccion focal Q k i el correspond iente codigo de alejamiento focal J ki = J/ en 
cada canal k. 

5 Este esquema resulta adecuado ya que, en general, las senales 
correspondientes a cada foco Megan a los elementos en momentos diferentes y 
su realizacion es sencilla. De este modo, la memoria de correccion focal tiene 
una anchura de palabra de b+1 bits, repartidos en dos campos: b bits para el 
codigo Q/c/y 1 bit para el J,. 

10 Una alternativa, que proporciona una mayor eficiencia, almacena en una 
memoria separada los codigos de alejamiento focal J„ cuyo contenido es 
compartido por un conjunto de canales. Para resolver el problema de no 
simultaneidad en la adquisicion de los focos, cada canal dispone una pequena 
cola donde se introducen los codigos J, obtenidos automaticamente de la 

15 memoria que los almacena. Cada canal utiliza el valor de J, en cabeza de la 
cola cada vez que alcanza un nuevo foco. Asf, una unica memoria central 
almacena un codigo J, por foco y unos pequenos FIFOs distribuidos facilitan la 
asignacion de estos codigos comunes a cada canal cuando lo necesita, como 
se expone con mas detalle en la descripcion del sistema. 

20 Por otra parte, puesto que la distancia entre focos puede modificarse 
dinamicamente, en la mayoria de las aplicaciones practicas basta una 
codificacion de la correccion focal con b=1 bit, lo que proporciona una alta 
eficiencia en el uso de memoria y una simplificacion de los circuitos de 
correccion. En particular, es inmediato realizar la distribucion uniforme de 0 6 1 

25 periodo de reloj maestro entre m muestras. 

Con lo expuesto, la Tecnica de Correccion Focal Progresiva, con su opcion 
Variable, proporciona un metodo que realiza la conformacion dinamica de 
haces en recepcion con una elevada resolucion temporal dentro de un gran 
rango de distancias. Adicionalmente, carece de redundancias en la adquisicion, 
30 no requiere procesos de interpolacion de muestras, permite mantener la tasa 
de muestreo mas baja compatible con el criterio de Nyquist y, por consiguiente, 
tiene un bajo consumo energetico. Su realizacion es compacta, debido tanto a 



WO 2005/073957 



PCT/ES2005/070008 



22 

un balance adecuado entre el calculo en tiempo real de las correcciones y el 
contenido de unas memorias precalculado en tiempo no real, como a que la 
codificacion de esta informacion es muy eficiente. 

Todo esto permite configurar un sistema versatil de conformacion de haces en 
5 recepcion con focalizacion dinamica, deflexion, apertura dinamica y 
apodizacion, cuyos detalles constructivos se exponen a continuacion. 

Descripcion del sistema. 

El sistema objeto de la presente invencion se basa en la tecnica de Correccion 
Focal Progresiva descrita anteriormente, incluyendo o no la opcion Variable, y 
10 realiza la focalizacion dinamica con deflexion de las senales recibidas por un 
conjunto de transductores o array, incluyendo funciones de apodizacion y de 
apertura dinamica. El sistema se describira con suficiente detalle para que 
pueda ser entendido y reproducido por personal con conocimientos de 
electronica digital. 

15 Con referenda a la Figura 2, un sistema de conformacion de haces en 
recepcion consta, en el caso mas general, de un array (10) compuesto de N 
transductores elementales (11) o, simplemente, elementos, que reciben 
senales del exterior, numerados 1, 2, ...k, N. Frecuentemente, estos 
elementos son, a su vez, emisores de senal para conformar un haz en emision, 

20 pasando a ser receptores una vez que esta finaliza. En la figura se han omitido 
los dispositivos necesarios para conformar el haz en emision. 

Las senales recibidas por los diferentes elementos son acondicionadas y 
amplificadas por los dispositivos AAS (12), que entregan las N senales 
analogicas ai, a 2 , a k , a N a otros tantos conversores analogico-digitales 

25 A/D (13). Los dispositivos AAS suelen incorporar, ademas de la funcion de 
amplificacion, otras que incluyen la compensacion de la atenuacion con la 
distancia mediante un control ganancia-tiempo, el filtrado y, en ocasiones, el 
desdoblamiento en las componentes en fase y cuadratura. En este ultimo caso, 
cada senal analogica a k esta formado por el par (i k , q k ) que representan, 

30 respectivamente, a las componentes en fase y en cuadratura. 
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Una vez acondicionadas, las senales analogicas son digitalizadas por los 
conversores A/D. En el caso de senales en banda base con pares (i k , q k ) 
pueden utilizarse dos conversores A/D por cada serial o multiplexarse un unico 
conversor A/D entre dos senales, las muestras pares para la componente i k y 
5 las impares para la q k o viceversa. 

Otra alternativa es la utilizacion de conversores AL, con los que las senales 
analogicas a k se representan mediante una secuencia continua de bits. En este 
caso la tecnica de Correccion Focal Progresiva y sus variantes resulta 
especialmente util al ajustar los instantes de muestreo a los de llegada de la 
10 serial a cada elemento con la resolucion de un reloj maestro cuyo periodo se 
adapta, ciclo a ciclo, al requerido por tales conversores AX para obtener una 
determinada relacion senal/ruido. 

Las tecnicas mencionadas son conocidas por los expertos en el tema y de 
aplicacion general. Las diferencias aparecen principalmente en el dispositivo 
15 CONF (14) o conformador de haz, cuya funcion es la de combinar las N 
secuencias numericas ei, e 2 , e k , e N proporcionadas por los conversores 
A/D para crear una unica secuencia de salida r que representa, en forma 
digital, a la serial recibida por el array focalizada, idealmente, en todos los 
puntos de una direccion determinada. 

20 Cada muestra entregada por los conversores A/D es una palabra de w bits 
siendo, tipicamente, 6 < w < 12 para el caso de conversores de muestreo 
instantaneo y 1 < w < 2 para los de tipo AZ, aunque el sistema puede operar 
con otras anchuras de palabra. 

Es, precisamente, la forma en que el dispositivo CONF realiza la funcion de 
25 combinar coherentemente las N entradas e k para generar la salida r lo que 
constituye el aspecto diferencial de la presente invencion. 

El dispositivo CONF (14) que es objeto de la presente patente, basa su 
funcionamiento en la tecnica de Correccion Focal Progresiva que constituye 
uno de sus aspectos claves, aportando diversas ventajas frente a otros tipos de 
30 realizaciones. 
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Asf, segun otros metodos resenados en la descripcion que se hizo del estado 
del arte, cada muestra de la secuencia numerica e k suministrada por k-esimo 
conversor A/D es sometida a un proceso de retardo como indica la Figura 3. 
Cada elemento se asocia a un canal de procesamiento individual (21). El 
5 contenido de una memoria de focalizacion dinamica (25) es utilizado por un 
dispositivo de control (24) para generar senales que producen un retardo 
grueso (22) y otro fino (23). El retardo grueso tiene, tfpicamente, la resolucion 
de un periodo del reloj de muestreo y el retardo fino una fraccion del mismo. 
Estos retardos se obtienen con diversos metodos que han sido descritos en la 
10 literatura especializada y son bien conocidos. 

El contenido de la memoria de focalizacion dinamica (25) es tal que, para cada 
muestra y para cada canal, los retardos introducidos alinean en el tiempo las 
salidas f l9 f 2 , fN que corresponden a los focos sucesivos F 19 F 2 , F N . Estas 
senales son sumadas en (20) para obtener la secuencia r focalizada 
15 dinamicamente. 

La mayor ventaja de estas tecnicas reside en que solo es necesario un reloj de 
muestreo, comun a todos los canales. Pero adolece de diversos inconvenientes 
derivados de la necesidad de realizar procesos de interpolacion para obtener 
una suficiente resolucion temporal en los retardos introducidos a las senales, 
20 con el consiguiente aumento de complejidad electronica y la posibilidad de 
introducir errores. Ademas, requieren una considerable cantidad de memoria 
para representar los retardos grueso y fino. 

El sistema objeto de la presente invencion opera de una forma diferente, 
muestreando directamente las senales en el instante en que Megan a cada uno 

25 de los elementos del array, procedentes de cada uno de los focos situados a lo 
largo de una direccion determinada. Al evitar la introduccion de retardos a las 
senales, se hacen innecesarios los procesos de interpolacion, eludiendo 
posibles errores y eliminando los circuitos electronicos correspondientes. 
Asimismo, solo se adquieren las muestras necesarias para realizar la 

30 composicion coherente, eliminando redundancias. 

En principio, el coste del metodo objeto de esta patente es que cada conversor 
A/D debe operar con un reloj de muestreo diferente, de frecuencia no uniforme. 
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Como se demuestra con la descripcion que sigue, la realizacion de esta 
caracterfstica es sencilla a partir de la metodologia enunciada como tecnica de 
Correccion Focal Progresiva, con la opcion Variable, pudiendo integrarse 
multiples canales en un unico dispositivo tipo FPGA estandar y comercialmente 
5 disponible. 

Con referenda de nuevo a la Figura 2, el metodo propuesto requiere que el 
dispositivo de conformacion de haz CONF (14), proporcione un reloj de 
muestreo independiente y de frecuencia no uniforme CK 1f CK 2 , CK k , CK N 
a cada conversor A/D (13), de forma que se adquieran las senales analogicas 

10 ai, a 2 , Bk, a N proporcionadas por los acondicionadores de senal AAS 
(12) en los instantes que corresponden a los de su llegada desde cada foco a 
cada elemento (11) del array (10). Las muestras obtenidas por los N 
conversores A/D (13) son procesadas para obtener su suma coherente, es 
decir, un valor de salida r que es la suma de las N muestras correspondientes a 

15 la senal recibida de un mismo foco por los N elementos del array. 

En la presente invencion, esta metodologia se aplica en forma modular, 
asociando un canal de procesamiento a cada elemento del array. La Figura 4 
muestra la estructura de uno de los canales (30) que procesa la senal recibida 
por el elemento k (11) del array. La estructura de todos los canales es identica 
20 a la mostrada. 

Como en otros casos, la senal es amplificada y acondicionada por el dispositivo 
AAS (12) por metodos conocidos, que incluyen dispositivos electronicos para 
amplificar, filtrar y manipular de forma analogica la senal. A la salida del AAS se 
obtiene la senal a k que es digitalizada por el conversor A/D (13). Las muestras 

25 son adquiridas en los instantes definidos por un generador de reloj (70) que 
produce una secuencia de pulsos CK k no uniformemente espaciados, cuya 
descripcion detallada se realiza posteriormente. El generador de reloj es 
habilitado por la senal HFD generada en el dispositivo HFOC (80) que marca 
los instantes en los que el canal comienza y termina la adquisicion de senal con 

30 focalizacion dinamica. 

La funcion del generador de reloj es la de garantizar que los instantes de 
muestreo determinados por CK k corresponden a los de la llegada de la senal 
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desde cada uno de los focos al elemento k, y opera segun los principios 
descritos anteriormente como Tecnica de Correccion Focal Progresiva, con su 
posible opcion Variable. 

Con esto, la secuencia de muestras e k obtenida a la salida del conversor A/D 
5 corresponde, precisamente, a los valores ordenados de la serial recibida por el 
elemento k procedente de cada uno de los focos. Estos valores son 
procesados por un bloque de apodizacion y apertura dinamica APD (60), que 
permiten mejorar las caracteristicas del haz y que produce la secuencia f k . Al 
existir una correspondencia 1:1 de las muestras e k a las f k , esta secuencia 
10 tambien esta ordenada con referenda a los instantes de llegada de la serial al 
elemento k desde cada foco. Los sucesivos valores de f k se almacenan 
transitoriamente en una cola o memoria FIFO que absorbe la diferencia en los 
instantes en que se producen muestras correspondientes a un mismo foco en 
diferentes canales. 

15 Elemento clave es, pues, el generador de reloj de muestreo (70), que produce 
la serial de digitalizacion CK k del k-esimo conversor A/D, de forma que las 
muestras obtenidas correspondan a los valores de la serial recibida por el 
elemento k desde cada uno de los focos. La Figura 5 muestra el diagrama de 
principio de una de las realizaciones preferidas del generador del reloj CK k , 

20 cuyo funcionamiento se basa en el metodo descrito anteriormente como 
Tecnica de Correccion Focal Progresiva, para el caso en que las correcciones 
focales se codifican con 6=1 bit. 

El circuito mostrado en la Figura 5 contiene solo los elementos basicos 
omitiendose otros dispositivos auxiliares, que son necesarios para el correcto 

25 funcionamiento, en favor de una mayor claridad de exposicion. Tiene un 
registro REG-A (71a) en el que se carga un valor relacionado con el parametro 
v o intervalo nominal entre focos en periodos del reloj maestro. Su salida V 
consta de multiples Ifneas que determinan el numero de etapas que componen 
el registro de desplazamiento de longitud programable SHR-A (72). Este 

30 dispositivo obtiene una salida p que es una version retrasada una cantidad de 
V periodos de reloj c de su entrada CK k . Este y otros dispositivos en el 
esquema mostrado operan smcronamente con el reloj maestro c de periodo T x . 
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La serial de entrada HFD habilita la operacion del registro de desplazamiento 
(72) y del flip-flop (77), de forma que estos dispositivos mantienen su estado 
actual independientemente del valor a su entrada cuando HFD=0. Cuando 
HFD=1, se habilita la focalizacion dinamica, interrumpiendose cuando HFD 
5 retorna a cero. Esta serial se genera en el dispositivo HFOC (80), Fig. 4. 

La senal p es, a su vez, retrasada 1 periodo de reloj maestro mediante el flip- 
flop FFA (73), produciendo la senal pr. El multiplexor MUX (74) proporciona una 
salida s que es, bien la senal p producida por el registro de desplazamiento o 
su version retrasada pr, en funcion del estado de la Imea de seleccion Q: 

10 s=p si Q=0; s=pr si Q=1 (35) 

El estado del selector Q lo determina el contenido de una memoria MEM (75) 
en la direccion F proporcionada por el contador CNT (76). En el caso actual, Q 
es una simple Imea (1 bit). 

La memoria se carga previamente, mediante circuitos no mostrados, con los 
15 codigos de correccion focal Q ki , donde k es el numero del canal actual y / es el 
numero de orden del foco, con /> 1. Los circuitos para la adquisicion de senal 
en el primer foco (#=0), que se describen posteriormente, determinan la 
activacion de la senal HFD. 

La salida del multiplexor (74) se registra en el flip-flop FFB (77), cuya salida 
20 CK k es el reloj de muestreo del conversor A/D del canal actual k (13), Figura 4, 
acciona al contador CNT (76) y es la entrada del registro de desplazamiento 
SHR-A (72). 

Para su funcionamiento, el registro de desplazamiento SHR-A (72) se inicializa 
a '0', el flip-flop FFA (73) a '0', el flip-flop FFB (77) a T y el contador (76) con la 
25 direccion de memoria donde se almacena el codigo de la correccion focal 
correspondiente al foco Fi (i=1). Los circuitos auxiliares para establecer estos 
estados iniciales se omiten en la figura. Asi, los valores iniciales de las 
diferentes senales son: CK k =1, p=0, pr=0 y s=0 (independientemente del 
estado de Q). 

30 A partir de este estado inicial, en el ciclo de reloj maestro c que sigue a la 
activacion de la senal HFD=1, se pone CK k =0 al encontrarse en este estado la 
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entrada s del flip-flop FFB (77) y el reloj habilitado por su entrada ce. Durante 
los V primeros periodos del reloj c se mantienen las senales p=0 y pr=0. En el 
ciclo siguiente se pone p=7, reflejando que el valor inicial CK k =1 se ha 
propagado por el registro de desplazamiento (72), pero todavfa se mantiene 
5 pr=0. En este mismo ciclo, si Q=0, es s=pr=0, pero si Q=1 , es s=p=1. Por 
consiguiente, en el siguiente ciclo, la salida CK k =1 solo si Q=7; tambien en 
este ciclo se pone pr=1 por lo que, si Q=0 se pondra la salida CK k =1 en el ciclo 
siguiente. Una vez que se hace CK k =1 el proceso se repite. 

Por consiguiente, el numero de ciclos n que transcurren desde el anterior valor 
10 CK k =1 y el proximo lo determina el valor Q de la siguiente forma: 

Si Q=0: n=V+1 

Si Q=7: n=V+2 (36) 

Acorde con la Ecuacion (20), cuando Q=0 el numero de ciclos entre focos es n 
= v y, cuando Q=7, debe ser n-v-1. Con esto: 

15 V=v-2 (37) 

que indica el valor inicial a cargar en el registro REG-A (71a) en funcion de la 
distancia nominal entre focos v expresada en periodos del reloj maestro. 

Cada vez que se produce un reloj de muestreo (CK k =1 ), se incrementa la 
direccion F proporcionada por el contador CNT (76), con lo que la memoria 
20 MEM (75) coloca en su salida el codigo correspondiente a la correccion focal a 
aplicar en el siguiente foco. 

Se observa que el conjunto FFA (73)-MUX (74)-FFB (77) configuran un registro 
de desplazamiento de longitud variable segun el codigo Q, que puede utilizarse 
como base para las realizaciones descritas a continuacion, aunque hay 

25 alternativas mas eficientes. Asf, la Figura 6 muestra un esquema de otra 
realizacion preferida que no limita el numero de bits con el que se representan 
los codigos de correccion focal Q. En este caso el conjunto FFA-MUX-FFB 
anterior se ha sustituido por un registro de desplazamiento SHRB (78) cuya 
longitud programable la determina el complemento del codigo de correccion 

30 focal Q obtenido mediante los inversores (75a) que producen las salidas NQ. 
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En este caso el intervalo entre los instantes de muestreo en dos focos 
consecutivos en el canal actual se determina por el valor multibit Q entregado 
por la memoria (75). Aquf el registro de desplazamiento SHRB (78) se inicializa 
aT. 

5 El funcionamiento es, por lo demas, analogo al caso anterior con la salvedad 
de que el numero de ciclos esta determinado por los valores V y NQ segun: 



Teniendo en cuenta que NQ=2 b -1-Q, el valor inicial a cargar en el registro REG- 
A (71a) en funcion del intervalo nominal entre focos en periodos del reloj 
1 0 maestro es: 



En una realizacion practica, los inversores (75a) son innecesarios si en la 
memoria de codigos de correccion focal (75), se almacenan los valores 
complementarios NQ en lugar de los de Q, manteniendose en la descripcion 
15 por coherencia en la notacion. Asimismo, en las realizaciones practicas de los 
circuitos mostrados en las Figuras 5 y 6 se incluye, ademas, logica auxiliar que 
evita la situacion en la que todas las salidas estan en estado '0' de la que, 
evidentemente, no saldrian los circuitos. 

La logica descrita permite focalizar cada una de las muestras adquiridas. 

20 Cuando se prefiera realizar una correccion focal periodica, compartiendo un 
mismo codigo de correccion focal Q cada m muestras, los circuitos anteriores 
se modifican como indica la Figura 7. El valor de m se carga en un registro 
RCM (71b) y el valor almacenado M es utilizado por una logica de control 
CTRL (90) junto con el valor Q del codigo de correccion focal entregado por la 

25 memoria (75) para generar las senales ng y cef. La primera controla la longitud 
del registro de desplazamiento SHRB (78) haciendola igual a mientras 
que cef habilita al contador CNT (76) que proporciona la direccion F a la 
memoria de correccion focal (75). 

El registro de desplazamiento SHRA (72) es de longitud programable en 
30 funcion del valor V programado en el registro REG-A (71a), con un numero de 



n=V+NQ 



(38) 



V=v -2 b +1 



(39) 
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etapas igual a V+1. El registro de desplazamiento SHRB (78) es tambien de 
longitud programable con un numero de etapas igual a ng+1. 

El numero de periodos de reloj maestro que median entre dos ciclos 
consecutivos CK k es igual al numero de etapas interpuestas entre la entrada de 
5 SHRA (72) y la salida de SHRB (78) que, acorde con lo anterior, es: 

n = V + ng + 2 (40) 

Programando en el registro REG-A (71a) un valor 

V=x>/m-3 (41) 

donde v/m representa el intervalo nominal entre muestras en periodos del reloj 
10 maestro, la Ecuacion (40) queda como: 

n =V /m + ng-l (42) 

La logica de control CTRL (90), mostrada en detalle en la Figura 7, controla el 
intervalo n mediante la serial ng, distribuyendo uniformemente el adelanto de 
periodos de reloj maestro Q entre las m muestras que se adquieren entre dos 
15 focos consecutivos. 

Para ello, el contenido CM de un contador modulo m, CMOD (91), que se 
actualiza con cada reloj CK k por su entrada de habilitacion ce, proporciona 
parte de la direccion DM de la memoria MDM (92). La otra parte de la direccion 
es el codigo de correccion focal Q. La direccion proporcionada a la memoria 
20 MDM es: 

DMU,QKJ-l>2 b +Q (43) 

donde j= CM+1 indica el numero de orden de la muestra actual. Programando, 
en tiempo no real, el contenido de las direcciones DM(j,Q) de la memoria MDM 
con los valores g(j,Q) que resultan de aplicar las Ecuaciones (22) a (24), la 
25 salida g representa el adelanto individual a aplicar a cada una de las muestras 
durante la adquisicion. El inversor (93) obtiene la serial ng, complementaria de 
la g, de modo que, conforme a la Ecuacion (42): 

Si g=0 ^ ng-1 =^> n - v/m 



WO 2005/073957 



PCT/ES2005/070008 



31 

Si g=1 =^> ng=0 =^> n = v/m - 1 (44) 

con lo que se produce un adelanto de 1 ciclo de reloj cuando g(j,Q) =1 y se 
mantiene el intervalo nominal v/m cuando g(j,Q) = 0, tal y como quedo descrito 
en el metodo de Correccion Focal Progresiva. 

5 Por su parte, cada vez que el contador CMOD (91) rebosa (esto es, conmuta 
del valor m-1 a 0), se produce un pulso cef que habilita el incremento del 
contador de focos CNT (76), proporcionando una nueva direccion F+1 a la 
memoria MEM (75). Asi, el proceso se repite con un nuevo valor de Q 
suministrado por la memoria MEM (75) con el contenido de la nueva direccion. 

10 Deben hacerse algunas observaciones adicionales. La primera es que, 
evidentemente, los inversores (93) se eliminan si en la memoria se programan 
los valores complementarios de los g(j,Q) dados por la Ecuacion (24). Otra es 
que la memoria MDM (92) es de pequeno tamano; por ejemplo, una memoria 
de 256x1 bits puede servir para valores m< 64 con b=2, o para m< 32 con b=3. 

15 A partir de este esquema se deriva la variante denominada como tecnica de 
Correccion Focal Progresiva Variable. Para ello, el registro RCM (71b) se 
sustituye por un contador, de forma que el numero de muestras m entre focos 
puede irse incrementando en el tiempo a partir de un valor mmimo inicial mi 
que, tfpicamente, puede ser igual a 1 . 

20 Tal y como fue descrito, segun la Ecuacion (18), la distancia minima Ro a partir 
de la cual el error absoluto del instante de muestreo se mantiene acotado en 
±T x /2 crece con el valor del intervalo entre focos v, expresado en periodos T x 
del reloj maestro. Al reducir la distancia entre focos en la region mas proxima al 
array mediante la programacion de un valor bajo para disminuye el valor de 

25 R 0 - Los instantes en que se incrementa el contador RCM (71b) vienen 
determinados por los instantes en que las senales proceden de focos situados 
a suficiente distancia como para que pueda adquirirse una muestra mas entre 
focos en todos los canales. El numero de orden de los focos en que esto ocurre 
se calcula en tiempo no real y se codifican en una tabla comun a todos los 

30 canales, a un subconjunto de canales o a un canal individual, estableciendo un 
balance entre eficiencia de uso de memoria y modularidad. 
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En una realizacion preferida y, en tiempo de adquisicion, los focos en los que 
se incrementa el contador RCM (71b) se determinan por el valor de un codigo 
de alejamiento focal J, expresado con 1 bit, con / = 1, 2, ... asociado a cada 
foco: el valor actual de m contenido en el contador RCM (71b) se incrementa si 
5 J, = 1 , manteniendo su valor actual si J, = 0 hasta el siguiente foco, donde se 
repite este proceso. El incremento de RCM solo se habilita, por consiguiente, 
cuando se produce el reloj de muestreo para adquirir la senal correspondiente 
a un foco. 

En una realizacion particular, la memoria de correccion focal (75) proporciona 
10 dos campos: uno, J ki = J, de un bit, indica si se incrementa el contador RCM 
(71b) en este foco y, el otro Qwde b bits, el adelanto a distribuir entre las m 
muestras actuales. En la figura se muestra la senal J=J, con una Ifnea de 
trazos. En este caso, preferentemente, se utilizan codigos de correccion focal 
Q/f/Con b=1 bit, puesto que puede controlarse de forma separada la distancia 
15 entre focos controlando el numero de muestras m intercaladas. En tal caso, la 
distribucion de los adelantos expresados por Q ki son de 0 6 1 ciclo de reloj 
maestro, realizandose en la muestra cuyo orden es la mitad del contenido de 
RCM (71b), aunque son posibles otras alternativas mas simples con poca 
incidencia en la calidad de la imagen. 

20 Sin embargo, para mejorar la eficiencia, en otra realizacion alternativa que 
tambien se muestra en la Figura 7, los codigos de alejamiento focal J, 
proceden de una memoria separada (MJ (59), Figura 10) cuya salida es 
compartida por multiples canales por la entrada J a una cola FIFO local (79). 
En tiempo de adquisicion, se realiza una lectura de un nuevo codigo de 

25 alejamiento focal J, en la memoria global MJ (59) de la Figura 10, cuando 
alguno de los canales a los que sirve detecta que su FIFO local (79) esta vacio. 
El codigo J, suministrado se escribe simultaneamente en todos los FIFOs 
locales (79) para ser utilizados por cada canal cuando adquieren senal de un 
nuevo foco, esto es, se activa su serial cef. En este momenta, el valor en 

30 cabeza de la cola Jc se utiliza para incrementar o no el valor de m contenido en 
el contador RCM (71b). 
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Debe destacarse que las colas FIFO locales (79) son de pequeno tamano, 
tfpicamente inferiores a 16x1 bit porque solo ha de almacenar un numero de 
codigos de 1 bit equivalente a la maxima diferencia de focos adquiridos en 
diferentes canales dentro de un mismo modulo. De este modo, con esta 
5 alternativa se consumen menos recursos que con la anterior en la que los 
codigos de alejamiento focal J, se almacenaban en la memoria de correccion 
focal (75) de cada uno de los canales. 

Con la Tecnica de Correccion Focal Progresiva y Variable la utilizacion de 
codigos de correccion focal se hace mas eficiente, distribuyendose entre un 
10 numero bajo de muestras al principio de la adquisicion (en el limite, 1 bit por 
muestra) y entre un numero mayor de muestras conforme la distancia de los 
focos al array aumenta (una fraccion decreciente de un bit por muestra). 



Por otra parte, es preciso determinar con precision el instante en el que se 
15 adquieren las muestras que corresponden al primer foco, de coordenadas (R 0 , 
Q), a partir del cual opera la correccion focal segun la metodologia descrita. 
Este instante esta determinado por la Ecuacion (19) para la geometria de la 
Figura 1 . Debe hacerse la observacion que el foco inicial puede situarse a una 
distancia mayor en funcion de las necesidades de la aplicacion. Tambien se 
20 tendra en cuenta la utilizacion o no del un numero variable de muestras entre 
focos para determinar este valor inicial. 

Esta funcion se realiza en el bloque HFOC (80) de la Figura 4, expuesto con 
mas detalle en la Figura 8. Una realizacion utiliza, simplemente, un contador 
accionado por el reloj maestro que se carga inicialmente con el valor N A (k): 



25 N A (k) = 



T 0 (k) 



(45) 

u 



donde T 0 (k) representa el intervalo entre el origen de tiempos hasta que la 
senal correspondiente al foco F 0 Mega al elemento k, determinado por la 
geometria y las velocidades de propagacion, o como expresa la Ecuacion (19) 
para la configuracion de la Figura 1 . 
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Las senales que corresponden al primer foco F 0 se adquieren, en el canal k, 
cuando rebosa este contador, habilitando, a partir de este momento, la 
actuacion de los circuitos de correccion focal progresiva. El error absoluto del 
instante de adquisicion de este primer foco es inferior a medio periodo de reloj 
5 maestro en todos los canales. Sin embargo, para mantener esta resolucion 
temporal, este contador debe ser accionado a la frecuencia del reloj maestro, 
normalmente elevada y muy superior a la de muestreo. 

Para evitar este inconveniente, una realizacion preferida, descompone el 
contador de ciclos N A en dos partes, como se muestra en la Figura 8, de forma 
10 que (omitiendo el mdice de canal k)\ 

N A =KN A1 +N A2 + 1 (46) 

con K, N A1 y N A2 enteros. El valor N A1 se carga en el contador CNT-A (82) 
mediante la serial Ida. Cuando se habilita su funcionamiento mediante la 
activacion de la serial cea, se acciona el contador con un reloj cb cu ya 
15 frecuencia es K veces inferior a la del reloj maestro, c, evitando las dificultades 
mencionadas. 

Tras N A1 ciclos cb, equivalentes a KN A1 periodos T x , el contador alcanza el 
valor terminal que alza la serial tca=1. En el siguiente ciclo de reloj cb se activa 
el flip-flop FF tipo RS (83), colocando la serial aj=1. Este valor se propaga por 

20 el registro de desplazamiento SHRG (84) cuya longitud la establece el valor 
programado para N A2 , accionado por el reloj maestro c. Transcurridos N A2 + 1 
periodos de reloj maestro, el valor aj=1 alcanza la salida del registro de 
desplazamiento, con lo que la salida pasa a HFD=1, habilitando a los circuitos 
de focalizacion dinamica una vez transcurridos N A periodos de reloj maestro, 

25 como indica la Ecuacion (46). De este modo, solo el registro de desplazamiento 
SHRG (84) debe operar a la elevada frecuencia del reloj maestro. 

Tal estado permanece estable hasta que se haya realizado un numero N F de 
correcciones focales, programado en el contador CNT-F (86). Cada correccion 
focal es indicada por la serial de entrada cef, que proporcion la logica de control 
30 CTRL (90) de la Figura 7 cuando se toman mutiples muestras entre focos, o es 
el reloj de muestreo CK k si solo se adquiere serial en los focos (Figuras 5 y 6). 
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Cuando CNT-F alcanza la cuenta terminal, se alza la senal tcf=1, provocando 
la puesta a cero del flip-flop FF (84) por su entrada R. Al realizarse la funcion 
AND (81), la senal HFD retorna a 0, finalizando el proceso de correccion focal 
progresiva. 

5 No obstante, la adquisicion de senales puede continuar con periodo de 
muestreo constante e igual al nominal v/m hasta completar el numero de 
muestras a adquirir, A/ s , programado en el contador CNT-S (88). El 
rebosamiento de este valor activa la senal FIN que termina por completo la 
adquisicion de senal en este canal. 

10 Es importante observar que, finalizado el proceso de correccion focal 
progresiva al hacerse HFD=0, la direccion F de la memoria de codigos de 
correccion focal (75) queda congelada en el ultimo valor alcanzado por el 
contador (76) de las figuras 5, 6 y 7, lo que posibilita iniciar automaticamente la 
siguiente adquisicion con los valores de correccion focal almacenados a 

15 continuacion, haciendo innecesario programar un nuevo valor inicial de F. De 
este modo es posible operar con diversas estrategias, como por ejemplo: 

a) Generacion de imagen sectorial, donde cada adquisicion se realiza con un 
angulo de deflexion diferente. Los codigos de correccion focal que 
corresponden a los angulos de deflexion sucesivos se almacenan 
20 secuencialmente en memoria, evitando tener que programar la direccion inicial 
F de la memoria de codigos de correccion focal antes de cada adquisicion. A 
diferencia de otras tecnicas, esta permite mantener el error de muestreo 
acotado y constante para todos los angulos de deflexion y a todas las 
distancias. 

25 b) Realizacion de barrido lineal u otras formas de adquisicion donde una 
subapertura activa compuesta por un subconjunto de los elementos del array, 
cambia en cada adquisicion, eliminando ciertos elementos e introduciendo 
otras. Analogamente, los codigos de correccion focal que corresponden a las 
configuraciones sucesivas se almacenan secuencialmente en memoria, 

30 evitando tener que cargar la direccion inicial F de la memoria de codigos de 
correccion focal antes de cada adquisicion. 
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c) Estrategias mixtas, en las que se adquiere senal de multiples direcciones 
utilizando los mismos codigos de correccion focal y diferentes valores N A (k) 
para diferentes direcciones con los que se adquiere el primer foco. En este 
caso, el valor de la direccion inicial debe restablecerse con anterioridad a la 
5 adquisicion. 

El dispositivo APD (60), Figura 4, realiza las funciones de apodizacion y 
apertura dinamica. La apodizacion, que consiste en multiplicar la senal recibida 
en cada canal por un numero real positivo A k comprendido entre 0 y 1 , mejora 
las caracterfsticas del haz conformado y el contraste en la imagen al reducir los 
10 lobulos laterales. La apertura dinamica permite mantener una resolucion lateral 
constante operando con un numero F# constante, que requiere que el tamano 
de la apertura activa aumente progresivamente en el tiempo. Ademas, la 
apertura dinamica facilita la aplicacion de la tecnica de Correccion Focal 
Progresiva en las regiones mas proximas al transductor. 

15 La Figura 9 muestra un diagrama de bloques del dispositivo APD (60). Para la 
funcion de apodizacion utiliza un multiplicador MULT (61) que obtiene el 
producto de cada uno de los valores de la secuencia e k por el coeficiente A k de 
apodizacion asignado al canal k. En una realizacion preferente, este coeficiente 
es un numero entero sin signo de a bits. Si los valores e k se expresan con w 

20 bits, de los a+w bits que componen el producto se mantienen los w bits mas 
significativos con redondeo en funcion del valor del primer bit ignorado y el 
signo de e k , incluyendose la logica correspondiente en el propio dispositivo 
MULT (61) para obtener la salida h k . La funcion realizada es: 



25 Debe hacerse la observacion de que la funcion de apodizacion puede 
habilitarse o no con sencillos circuitos de control no mostrados en la figura 9. 
Basicamente, la salida h k es la entrada e k si la funcion de apodizacion esta 
inhabilitada, respondiendo a la Ecuacion (47) en otro caso, funcion que se 
realiza con un simple multiplexor gobernado por la senal de habilitacion de la 

30 apodizacion. 




JU 



(47) 
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Para la funcion de apertura dinamica, se dispone un contador binario CNZ (62) 
que se carga inicialmente con el valor MZ(k) o numero de muestras nulas 
expresado por la Ecuacion (25) para el canal k. Este contador se acciona con 
un reloj cuya frecuencia es la nominal de muestreo, cs. Su operacion se habilita 
5 con la activacion de la senal HFD=1 que inicia el proceso de focalizacion 
dinamica y es proporcionada por el dispositivo HFOC (80), Figuras 4 y 8. Por 
su parte, la senal HFD=1 actua sobre la senal S del flip-flop FF (64) poniendo la 
senal z=7 que selecciona la entrada nula (63) al multiplexor MUX (65). 

En este estado, se inicia la adquisicion y procesamiento de senales en todos 
10 los canales. En aquellos en los que el contador CNZ (62) no haya alcanzado la 
cuenta terminal (senal fc), las muestras de salida en f k seran de valor nulo, por 
lo que el sistema opera como si el elemento asociado no estuviera active Tras 
MZ ciclos del reloj de muestreo, durante los que se han producido otras tantas 
muestras nulas sobre f k , se alza la senal tc=1 que pone a cero el flip-flop FF 
15 (64) y habilita la salida de los valores h k obtenidos por el multiplicador (61). De 
este modo se proporciona un numero programable de muestras nulas al 
comienzo de cada adquisicion en cada canal, lo que permite programar las 
caracteristicas de la apertura dinamica. La secuencia f k es: 

f k = {0, 0, ...,0, Ai<e k (NZ+1), A^NZ+2), ...} (45) 

20 donde e k (i) es la muestra / de la secuencia e k y que, por consiguiente, se 
corresponde con la muestra / de cualquier otra secuencia. De este modo se 
mantiene el orden para realizar posteriormente la composicion coherente de 
todas las senales, sean o no nulas. 

25 Otro de los aspectos importantes de la presente invencion es su capacidad de 
realizacion modular, cuyo concepto se muestra en la Figura 10. Una realizacion 
preferente agrupa en un modulo (50) todas las funciones de procesamiento de 
senal para un subconjunto de PS elementos del array, asf como otras 
funciones de control global. En particular, un modulo incluye la memoria MJ 

30 (59) que proporciona los codigos de alejamiento focal J, junto con los circuitos 
necesarios para su lectura y carga inicial. Ademas, cada modulo (50) contiene 
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S submodulos (40), donde cada submodulo procesa las senales de P 
elementos, y un combinador (55) que permite enlazar multiples modulos en 
cadena. 

Internamente, los S submodulos (40) que conforman un modulo, se encadenan 
5 para realizar la composicion coherente de los P S canales en cascada a traves 
de los FIFOs (52). Esta estructura facilita la operacion en paralelo de todos los 
canales de procesamiento digital, donde los FIFOs compensan 
automaticamente las diferencias en los instantes de adquisicion de serial en los 
diferentes canales. 

10 Para ello, cada uno de los elementos (11) del array (10) proporciona la senal 
recibida a uno de los P procesadores elementales, PE, (30) integrados en un 
submodulo (40). Cada procesador elemental PE realiza el acondicionamiento 
de senal y digitalizacion en el instante de llegada de la senal al elemento que le 
corresponde, siguiendo la metodologia descrita anteriormente, e incluye los 

15 dispositivos que fueron mostrados con detalle en las Figuras 4 a 8. 

Cada PE proporciona a su salida una secuencia ordenada de valores 
correspondientes a las senales recibidas por el elemento asociado desde cada 
uno de los focos o muestras intercaladas. Estas salidas se suman (41) para 
realizar la composicion coherente de las P senales que corresponden a un 
20 submodulo. 

Por su parte, el dispositivo (42) suma los resultados obtenidos localmente en el 
submodulo con los que le haya colocado el precedente en su cola FIFO (52). 
De este modo, a la salida de cada submodulo se obtiene una secuencia 
ordenada que corresponde a las senales focalizadas de sus P elementos mas 
25 los de los submodulos que le preceden en la cadena, en el orden de llegada de 
senales desde cada foco. El primer submodulo de la cadena recibe un valor '0' 
(51) en sustitucion del FIFO (52), considerando que no existe un submodulo 
precedente. 

A su vez, el combinador (55) suma en (58) los resultados obtenidos por este 
30 modulo y que se almacenan en el FIFO (56) con los que el modulo anterior ha 
colocado en la cola FIFO (57). Como en ambos casos los resultados son 
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secuencias ordenadas de valores correspond ientes a focos sucesivos, la 
secuencia de salida del sumador refleja este mismo orden. El encadenamiento 
de L modulos (50), cada uno de ellos con S submodulos (40) con P 
procesadores elementales en cada submodulo proporciona un sistema de 
5 focalizacion dinamica con la tecnica de Correccion Focal Distribuida para 
N=LSP elementos. 

Es, por tanto, una estructura que presenta diversas ventajas: 

• Es modular, facilitando el desarrollo de circuitos y mayor flexibilidad de 
aplicacion. 

10 • Procesa en paralelo, lo que permite alcanzar elevadas tasas de calculo. 

• Tiene una capacidad ilimitada de expansion desde el punto de vista logico. 

La Figura 11 muestra una de las posibles realizaciones del sumador (41) 
integrado en un submodulo (40). Este sumador recibe las P secuencias ei a e P 
correspond ientes a los valores de las muestras adquiridas por los P 
15 conversores A/D pertenecientes al sub-modulo. En el caso mostrado en la 
figura es P=8, siendo la estructura generalizable a cualquier otro valor, 
preferiblemente potencia de 2. 

Las muestras que corresponden a cada foco en las secuencias a e P se 
obtienen en forma ordenada, esto es, la primera muestra obtenida por los P 

20 conversores A/D corresponde al primer foco, la siguiente al segundo foco y asi 
sucesivamente. Este razonamiento es extensivo aplicado a las muestras 
obtenidas entre focos, con o sin la opcion de Correccion Focal Progresiva 
Variable. Sin embargo, los instantes de muestreo no son simultaneos en todos 
los canales, sino que es posible que un canal adquiera varias muestras 

25 sucesivas mientras que otros esperan a que las senales correspondientes 
lleguen al elemento asociado. 

Este desfase se corrige automaticamente con las memorias FIFOs distribuidas 
(43), (45) y (47) junto con los sumadores (44), (46) y (48). El FIFO (49) realiza 
la misma funcion pero con referenda a las secuencias de muestras de dos 
30 submodulos. 
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Los FIFOs se asocian en parejas, cuyas salidas se suman cuando ambos 
disponen de, al menos, una muestra almacenada. Esta situacion se detecta 
con un sencillo circuito de control que ejecuta una lectura simultanea en ambos 
FIFOs cuando no estan vacios y escribe el resultado obtenido por el sumador 
5 correspondiente en el FIFO conectado a su salida. Cada lectura libera una 
posicion de memoria en los FIFOs, que puede ser ocupada con la escritura de 
una nueva muestra. Cada escritura ocupa una nueva posicion de memoria en 
el FIFO correspondiente. Las lecturas en los FIFOs asociados son simultaneas, 
pero no asi las escrituras que dependen de los instantes de muestreo en los 
10 respectivos conversores A/D, a su vez derivados del instante de llegada de la 
senal desde cada foco al elemento. 

La capacidad H de cada FIFO (43) debe ser suficiente para almacenar las 
muestras adquiridas hasta que se disponga de una muestra en el canal 
asociado, en cuyo momento opera el sumador (44), liberando la posicion 
15 ocupada por las muestras correspondientes y transfiriendo el resultado de la 
suma al FIFO (45). 

Analogamente, los FIFOs (45) deben tener una capacidad 2-H, los FIFOs (47) 
una capacidad 4 H, y asf sucesivamente. De este modo, la cantidad de 
memoria distribuida en forma de FIFOs requerida por un submodulo con P 
20 entradas es: 

HP + 2HPI2 + 4H*P/4 + ...= H P log 2 P (49) 

En el caso en que se utilizara un unico FIFO por cada canal, su capacidad 
deberia ser de H(P-1) posiciones, con lo que el total de memoria requerida 
seria HP(P-1). Asi, por ejemplo, si H=16y P=8, segun el metodo preferido 
25 cuya estructura se muestra en la Figura 6 se requieren 384 posiciones de 
memoria distribuida, mientras que en el caso de utilizar un unico FIFO por cada 
canal serian necesarias 896 posiciones. 

El FIFO (49), por su parte, almacena transitoriamente los resultados obtenidos 
por el submodulo actual hasta que el sumador (42) de la Figura 10 pueda 
30 liberar una posicion cuando el FIFO (52) disponga de un dato correspondiente, 
entregando en tal caso el resultado al submodulo siguiente. 
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Breve description del contenido de las figuras. 

La Figura 1 muestra la geometria utilizada como ejemplo, con un array 
compuesto por N elementos numerados 1, 2, 3, ... k, AM y A/, sobre el eje x. 
Los focos se situan a intervalos AR a partir de una distancia minima R 0 . Se 
5 muestra el segmento entre un foco situado en coordenadas polares (R, Q) y el 
elemento k del array situado en coordenadas cartesianas (x k , 0), cuya longitud 
es L k (R,Q). 

La Figura 2 muestra la estructura general de un sistema de conformacion digital 
de haces en recepcion generico. Se compone de un array (10) con N 
10 elementos (11) y un canal de procesamiento de serial por cada elemento. Cada 
canal contiene circuitos de amplificacion y acondicionamiento de serial AAS 
(12) y un digitalizador A/D (13). Las salidas de los N digitalizadores son 
procesadas por el conformador CONF (14). 

La Figura 3 muestra la estructura utilizada frecuentemente para realizar la 
15 composicion coherente de senales segun metodos conocidos. Consta de un 
sistema de retardo grueso (22) y otro fino (23), cuya operacion se gestiona 
desde un controlador local (24) que utiliza el contenido de una memoria de 
focalizacion dinamica (25). Las N salidas retardadas f 1t f 2 , se suman en 

(20) para obtener la secuencia resultante r. 

20 La Figura 4 muestra la estructura de un canal de procesamiento digital 
preferentemente utilizado por el sistema objeto de esta patente. Contiene el 
elemento de amplificacion y acondicionamiento de serial AAS (12), el 
digitalizador A/D (13), una memoria FIFO (43) que almacena transitoriamente 
las muestras adquiridas, un generador de reloj de muestreo (70) que constituye 

25 un elemento fundamental de la presente invencion y un modulo de habilitacion 
de la focalizacion dinamica y control, HFOC (80). 

La Figura 5 muestra la estructura de principio para realizar el generador del 
reloj de muestreo (70) con codigos de correccion focal Q expresados con 1 bit 
en la memoria MEM (75). El contador CNT (76) establece la direccion de 
30 acceso a esta memoria. El reloj de muestreo CK k , propio del canal /c, se genera 
mediante logica distribuida en el registro REG-A (71a), el registro de 
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desplazamiento de longitud programable SHR-A (72), el flip-flop FFA (73), el 
multiplexor MUX (74) y el flip-flop FFB (77). 

La Figura 6 muestra la estructura de principio para realizar el generador del 
reloj de muestreo (70) con codigos de correccion focal Q expresados con 
5 multiples bits en la memoria MEM (75). El contador CNT (76) establece la 
direccion de acceso a esta memoria. El reloj de muestreo CK k , propio del canal 
k, se genera mediante logica distribuida en el registro REG-A (71a), el registro 
de desplazamiento de longitud programable SHRA (72), el registro de 
desplazamiento de longitud programable SHRB (78) y el conjunto de inversores 
10 (75a). 

La Figura 7 muestra la estructura de principio para realizar el generador del 
reloj de muestreo (70) con codigos de correccion focal Q expresados con uno o 
mas bits en la memoria MEM (75). El contador CNT (76) establece la direccion 
de acceso a esta memoria. En este caso puede insertarse un numero fijo o 

15 variable de muestras entre focos, segun el contenido del registro-contador 
RCM (71b) controlando la variable ng que establece la longitud del registro de 
desplazamiento SHRB (78). La variable ng, de 1 bit, es proporcionada por el 
control CTRL (90) mostrado con mas detalle en la parte inferior de la figura. 
Esta unidad consta de un contador de modulo programable CMOD (91), cuyo 

20 contenido y el valor del codigo de correccion focal Q, forman la direccion a una 
memoria local MDM (92). La salida g de esta memoria es invertida por (93) 
para obtener la senal ng. En este caso, existe la opcion de proporcionar el 
codigo de alejamiento focal J desde la misma memoria MEM de codigos de 
correccion focal (75), como se muestra con trazos de puntos o, 

25 preferiblemente, desde el exterior y a traves de una memoria FIFO (79). Por lo 
demas, el reloj de muestreo CK k del canal k se produce mediante la logica 
compuesta por el registro REG-A (71a), el registro de desplazamiento SHRA 
(72) y el registro de desplazamiento SHRB (78). La unidad de control CTRL 
proporciona, ademas, la salida cef utilizada para habilitar el conteo de focos en 

30 la unidad HFOC (80) mostrada en las Figuras 4 y 8. 

La Figura 8 muestra esquematicamente el principio en que se basa la unidad 
HFOC (80) que realiza el control general de la adquisicion. Consta de un 
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contador CNT-A (82) que se carga mediante Ida con el valor N A i, se habilita 
con la senal cea y se acciona con el reloj de baja frecuencia cb. La cuenta 
terminal tea actua sobre la entrada S del flip-flop FF (83). Su salida ajse 
retrasa una cantidad de ciclos de reloj maestro c expresada por N A 2 mediante 
5 el registro de desplazamiento SHRG (84). La puerta AND (81) produce un nivel 
alto en la salida HFD cuando sus dos entradas estan en alto. El flip-flop (83) se 
repone a 0 cuando se alza su entrada R al activarse la senal tcf. 

La senal tcf es la cuenta terminal del contador CNT-F (86) que puede 
inicializarse con un numero de focos A/Fdeterminado al activarse la senal Idf, se 
10 habilita con la senal cef y se acciona con el reloj de baja frecuencia cb. Por su 
parte, el contador CNT-S (88) cuenta las muestras adquiridas, activando su 
salida FIN cuando se alcanza el numero programado Afe. Esta programacion se 
realiza activando la senal Ids, el contador se habilita con ces y se acciona con 
el reloj de baja frecuencia cb. 

15 La Figura 9 muestra esquematicamente el principio en que se basa la unidad 
APD (60) que realiza las funciones de apodizacion y apertura dinamica. Consta 
de un multiplicador (61) que obtiene el producto de la entrada e k por el del 
coeficiente A k y entrega su resultado h k a un multiplexor MUX (65). Asimismo, 
dispone de un contador CNZ accionado por un reloj cs y habilitado por la senal 

20 HFD cuyo contenido inicial es M z , esto es, el numero de muestras nulas a 
aportar por este canal para realizar la funcion de apertura dinamica. Al alzarse 
la senal HFD se pone en alto la senal z de salida del flip-flop FF (64), 
seleccionandose la salida '0' entregada por la unidad (63) por el multiplexor 
MUX (65), que la proporciona como salida f k . Cuando se alcanza la cuenta 

25 terminal tc se activa la entrada R del flip-flop FF (64), pasando a nivel bajo su 
salida z con lo que el multiplexor MUX (65) selecciona la salida h k 6e\ 
multiplicador MULT (61), entregandola por la salida f k . 

La Figura 10 muestra esquematicamente la arquitectura modular del sistema 
objeto de esta patente. Un modulo (50) contiene multiples submodulos (40), 
30 cada uno de ellos con varios procesadores elementales PE (30) cuyas entradas 
son las senales recibidas por los elementos (11) del array (10). Las salidas de 
los procesadores elementales son sumadas en (41) y la salida de este sumador 
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se utiliza como uno de los sumandos de una cadena de sumadores (42) y 
FIFOs (52). El primer submodulo sustituye este FIFO por una entrada nula (51). 
Los resultados proporcionados por el sumador (42) del ultimo submodulo se 
almacenan transitoriamente en un FIFO (56) para ser, a su vez, sumados en 
5 (58) con los resultados de un modulo precedente que han sido almacenados en 
el FIFO (57) en el sumador (58). Adicionalmente, el modulo contiene una 
memoria MJ (59) que proporciona los codigos de alejamiento focal J a los 
diferentes procesadores elementales PE (30). 

La Figura 11 muestra esquematicamente el principio en que se basa una 
10 realizacion preferida del sumador (41), con un ejemplo para 8 canales, donde 
multiples FIFOs (43) almacenan transitoriamente las muestras obtenidas en 
cada procesador elemental e k . Las salidas de los FIFOs (43) son, a su vez, 
sumadas por pares en los sumadores (44), almacenandose el resultado en los 
FIFOs (45). Las salidas de estos FIFOs son sumadas en (46) y el resultado se 
15 almacena en los FIFOs (47). Analogamente, las salidas de estos FIFOs son 
sumadas en (48), almacenandose los resultados finales en el FIFO (49), cuya 
salida r representa la suma de las muestras ei a e 8 entregadas por los 8 
procesadores elementales que constituyen un submodulo en este ejemplo. 

La Figura 12 muestra un esquema de un ejemplo realizado de un modulo (200) 
20 del sistema objeto de esta patente. Dispone de conectores (120a) y (120b) por 
los que se introducen las senales analogicas que son distribuidas por los buses 
(121a) y (121b) a los conversores A/D dobles (101). Las muestras 
proporcionadas por estos conversores son introducidas por (106) en una FPGA 

(102) que realiza las funciones de procesamiento digital correspondientes a un 
25 submodulo (100), entre las que destacan la generacion de los relojes de 

muestreo (105) y el suministro de resultados al siguiente submodulo por el bus 

(103) . El ultimo submodulo en la cadena suministra los resultados obtenidos a 
una FPGA (111) por el bus (104) para ser combinados con los resultados 
procedentes de otro modulo por el bus (113), suministrandose los datos 

30 resultantes por los buses (112) o (114). Estos buses estan presentes en el 
conector (110). A su vez, la FPGA (111) entrega senales de control (115a) y 
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(115b) a los circuitos de amplificacion y acondicionamiento analogico a traves 
de los conectores (130a) y (130b), respectivamente. 

La Figura 13 muestra las medidas experimentales realizadas sobre el modulo 
descrito en la Figura 12, mostrando el error absoluto en ns para cada uno de 
5 los 32 canales. 

EJEMPLO DE REALIZACI6N DE LA INVENCI6N: 

En todos los casos descritos se obtiene el reloj de muestreo de cada canal 
utilizando muy pocos recursos hardware, fundamentalmente registros, 
10 contadores y una memoria para almacenar los codigos de correccion focal, 
realizandose la adquisicion de las senales con una alta resolucion temporal. 

Debido a la alta eficiencia de la metodologia descrita en cuanto a utilizacion de 
recursos hardware y, a diferencia de otras realizaciones, las tecnicas CFP y 
CFPV permiten integrar multiples canales de procesamiento en un unico 
15 circuito VLSI, junto con otras funciones descritas, incluso en dispositivos de uso 
general como son las FPGAs. 

Asf, en una realizacion demostrativa de esta invencion, un modulo contiene 4 
submodulos, cada uno de ellos con 8 procesadores elementales. Este pequeno 
sistema puede operar con una apertura activa de hasta 32 elementos en 
20 barrido lineal o sectorial con apertura dispersa y arrays de 128 elementos. Con 
4 modulos se configura un sistema que opera con aperturas densas o 
aperturas activas de hasta 128 elementos. El sistema es escalable, anadiendo 
32 elementos activos por cada nuevo modulo instalado. 

La Figura (12) muestra graficamente la configuracion del modulo, realizado en 
25 una tarjeta de circuito impreso de 220x100 mm. En la figura se han omitido los 
circuitos correspond ientes al preprocesamiento analogico de senales, con 
posible excitacion de los transductores, asf como los de posprocesamiento que, 
en esta realizacion particular, se basan en un computador integrado al que se 
conecta una Unidad de Control e Interfaz (UCI) con el sistema descrito. 
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El modulo (200) recibe las senales analogicas de 32 canales proporcionadas 
por los circuitos de amplificacion y acondicionamiento por los conectores (120a) 
y (120b). Contiene cuatro submodulos (100) y un combinador realizado en la 
FPGA (111). El combinador recibe senales procesadas por un modulo anterior 
5 en la cadena por el bus (113) y envfa los resultados al que le sucede por el bus 
(114). Los parametros se programan desde la UCI a traves del bus de control 
(112) por el que, tambien, el ultimo modulo de la cadena entrega los resultados 
finales. Asimismo, la FPGA (111) entrega un conjunto de senales digitales por 
los conectores (130a) y (130b) que son utilizadas por la unidad de 
10 procesamiento analogico. Por lo demas, esta estructura refleja fielmente la 
arquitectura presentada en la Figura 10. 

Las senales analogicas introducidas por los conectores (120a) y (120b) se 
distribuyen a los submodulos mediante los buses (121a) y (121b), 
respectivamente, donde son digitalizadas por los conversores A/D (101). La 
15 salidas digitales de los conversores A/D se introducen en las FPGAs (102) 
mediante los buses (106). Los relojes de muestreo de dichos conversores (105) 
se generan en las FPGAs (102) con la tecnica de Correccion Focal Progresiva 
descrita en esta memoria. 

Por su parte, los submodulos se conectan entre sf mediante los buses 
20 bidireccionales (103). Junto al bus (104) que conecta el ultimo submodulo al 
combinador (111), estos buses se utilizan en un sentido para la programacion 
de parametros desde el de control (112) y en el sentido opuesto para 
suministrar resultados a los submodulos posteriores en la cadena o al 
combinador (111). 

25 El combinador realiza la suma coherente de los resultados obtenidos en el 
modulo actual a traves del bus (104) con los que recibe de un modulo anterior 
por el bus (113), suministrando los resultados de este proceso por el bus (114) 
al modulo siguiente o por el bus (112) a la Unidad de Control e Interfaz. Una u 
otra opcion son programables. 

30 En esta realizacion se han utilizado 16 conversores dobles de 10 bits, tipo 
AD9218 de Analog Devices para los dispositivos (101), 4 FPGAs tipo XC2S200 
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de Xilinx Inc. para los dispositivos (102) y una XC2S150 de Xilinx Inc. para el 
dispositivo (111). 

El reloj del sistema es de 40 MHz, velocidad a la que operan smcronamente las 
comunicaciones entre submodulos, modules y con la UCI. Internamente, parte 
5 de los circuitos tambien operan a esta frecuencia o una fraccion que determina 
la de muestreo. Sin embargo, en las FPGAs (102) se multiplica por 4 la 
frecuencia del reloj del sistema para obtener el reloj maestro de 160 MHz. La 
resolucion temporal es, pues, T x = 6.25 ns. y el error de muestreo esta 
teoricamente acotado a I £ I < 3. 125 ns. Con esto, pueden procesarse senales 
10 analogicas de hasta 10 MHz con jli > 16 y hasta 20 MHz con [i> 8, siendo jli la 
relacion entre el periodo de la senal y la resolucion temporal que afecta al 
rango dinamico y contraste de las imagenes. 

Cada FPGA (102) contiene los circuitos necesarios para realizar todo el 
procesamiento digital de 8 canales, incluyendo: 

15 • Apodizacion, con coeficientes 0 < A k < 255/256 o deshabilitada (A k = 1, V k). 

• Apertura dinamica, con 0 < MZ k < 2 12 . 

• Deflexion del haz con foco fijo, con 0 < NA k < 2 14 . 

• Focalizacion dinamica segun la tecnica de Correccion Focal Progresiva 
expuesta, con 4096 codigos de correccion focal de 1 bit por cada canal. 

20 • Programacion del numero de muestras a adquirir, con 2 < Ns < 2 17 . 

• Programacion del numero de muestras entre focos, con 4 < m < 64. 

• Programacion del numero de focos por adquisicion, con 0 < N F < 2 12 . 

• Programacion de la frecuencia de muestreo de salida, con 2.5 < fs < 40 
MHz. 

25 • Circuitos de gestion de ahorro energetico y habilitacion de funciones. 

• Circuitos de programacion de parametros y procesamiento de resultados. 

• Circuitos para la verificacion del modulo, integrando un generador digital de 
serial. 
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Por su parte, el combinador contenido en la FPGA (111) tiene circuitos para 
realizar: 

• La gestion del modulo, transfiriendo valores de parametros a los 
5 submodulos por (104) y (103). 

• La recopilacion de resultados procedentes del modulo anterior por (113). 

• La composicion de estos resultados con los del propio modulo, recibidos por 
(104). 

• El suministro de resultados combinados al modulo siguiente por (114). 
10 • El suministro de resultados finales a la UCI a traves del bus (112). 

• Circuitos para compensar el offset de las senales. 

• Circuitos para controlar la unidad de pre-procesamiento analogico. 

La Figura (13) muestra medidas experimentales relativas al error absoluto en 
15 ns. (ordenadas) del instante de muestreo en cada uno de los 32 canales (en 
abcisas) para un foco situado a una distancia de 50 mm., con un angulo de 
deflexion de 15°, velocidad de propagacion c= 5900 nrvs' 1 y array lineal de 5 
MHz con elementos uniformemente distribuidos a intervalos 6=0.8 mm. 

El maximo error teorico es de 3.125 ns., que se verifica en 28 de los 32 
20 canales. El mayor error en los 4 canales que superan el valor teorico se debe a 
diferencias en el trazado de pistas en el circuito impreso y en las FPGAs de 
este modelo experimental y, en cualquier caso, no es significative El error 
cuadratico medio, mas determinante de la calidad de la imagen, es de 1.9 ns. 

Es evidente que pueden introducirse mejoras a la realizacion descrita puesto 
25 que existen dispositivos con mayor capacidad de integracion. Asimismo, 
pueden realizarse otras modificaciones que afectan al tipo de conversion A/D, 
distribucion de dispositivos, etc., ya que el modelo descrito no es mas que un 
ejemplo demostrativo, sin que tales modificaciones supongan un cambio 
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sustancial a los fundamentos de la metodologia y sistema que han sido 
descritos en la presente memoria. 

Para la aplicacion del metodo, considerese un caso de imagen ultrasonica 
sectorial en pulso-eco, con un transductor lineal de N elementos centrado en el 
5 origen, con un tamano de apertura D. Se pretende cubrir con la imagen un 
sector de -45° < 9 < 45°. 

Si fR es la frecuencia de la senal recibida, para verificar con suficiente margen 
el criterio de Nyquist, se elige una frecuencia de muestreo nominal fs = 4 f R . 
Ademas, para mantener un elevado rango dinamico y contraste en la imagen, 
10 se elige |u = fx/fe = 16, siendo f x la frecuencia del reloj maestro. 

Considerando el trayecto en ida y vuelta del pulso ultrasonico, el intervalo entre 
muestras es: 

AR = cT s /2 (50) 
Sustituyendo el valor AR dado por la Ecuacion 1 resulta, con m = 1: 

15 V=|f = 4 (51) 

Cuando se insertan m muestras entre focos, el valor de v = 4 m. Aplicando 
ahora la Ecuacion 18, con \x k \ < D/2, donde D es el tamano de la apertura, la 
distancia mmima para aplicar la Tecnica de Correccion Focal Progresiva con 
apertura completa es: 



20 



> — 



D 2 



mv -2a 

COS0 



2^amv - a 2 



+ ^|sen6| (52) 



donde a=2 b -1, siendo b el numero de bits con los que se expresan los codigos 
de correccion focal, con la condicion mv > 2a. 

La Tabla II muestra los resultados de la Ecuacion (52) para b-1, 2 y 3 en 
funcion del numero de muestras (fijo) entre focos. 

25 Debe observarse que el valor Ro/D representa el mmimo valor del numero F# 
con el que se puede operar utilizando la funcion de apertura dinamica. En 
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general, se acepta como valido un valor F# min < 1 , que se obtiene con la 
mayoria de las combinaciones. 

Tabla II 

Relacion F# min = Ro/D para aplicar la TCFP con apertura completa (-45°< 9 < 
5 45°) 



m 


b = 1 


b = 2 


b = 3 


1 


0.6 






2 


0.8 


0.4 




4 


1 


0.6 


0.4 


8 


1.4 


0.9 


0.6 
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REIVINDICACIONES 

Esta patente descubre el metodo y sistema para conformacion de haces en 
recepcion con una nueva tecnica denominada de Correccion Focal Progresiva, 
con la opcion Variable. Segun dichos metodo y sistema, se realiza la 
5 composicion coherente de senales procedentes de fuentes situadas a lo largo 
de una direccion seleccionada, que Megan a un conjunto de A/transductores 
receptores o elementos distribuidos espacialmente, formando un array. Un 
conformador de haces en recepcion realiza una composicion coherente de las 
senales recibidas que proceden de fuentes puntuales situadas a lo largo de 

10 dicha direccion. La composicion coherente de las senales originadas en cada 
punto y recibidas por cada uno de los N elementos del array proporciona una 
unica senal de salida que representa a la que se obtendria con un transductor 
receptor unico, del mismo tamano que el array, focalizado individualmente en 
todos y cada uno de los focos, tecnica que se conoce en el ambito como 

15 focalizacion dinamica con deflexion. 

Las senales recibidas son amplificadas y acondicionadas analogicamente antes 
de ser digitalizadas. La tecnica de correccion focal progresiva, con o sin la 
opcion variable, y el sistema que la realiza descritos en esta patente, se aplican 
a las senales digitalizadas por un conversor analogico-digital. Se denomina 
20 canal al conjunto de dispositivos que procesan la senal recibida por un 
elemento del array hasta su combinacion con las correspondientes a otros 
elementos. 

A diferencia de otros metodos que forman parte del estado del arte, el metodo y 
sistema objeto de esta patente no se basa en retrasar las senales recibidas por 

25 cada uno de los elementos del array, ni utiliza procesos de interpolacion. Antes 
bien, se fundamenta en muestrear las senales recibidas en los instantes en que 
Megan a cada uno de los elementos del array desde cada uno de los focos. El 
metodo garantiza que la composicion de muestras es coherente, esto es, solo 
se combinan las que corresponden al mismo foco y, ademas, no redundante, 

30 adquiriendo unicamente las que van a utilizarse en el proceso de composicion. 
Requiere un reloj de muestreo de las senales no uniforme y diferente para cada 
uno de los canales. 
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Por otra parte y, a diferencia de otras tecnicas en las que cada reloj de 
muestreo queda establecido por una secuencia de '1' y '0' en una memoria 
individual para cada canal o se calcula en tiempo real mediante circuitos 
digitales especificos a partir de los valores de una diversidad de parametros 
5 programados durante la adquisicion, con el metodo y sistema objeto de esta 
patente los instantes de muestreo se determinan interpretando el contenido de 
tablas asociadas a cada canal, que se calculan y almacenan en memorias 
antes del inicio de la adquisicion de serial, haciendo innecesaria la carga de 
multiples parametros en tiempo de adquisicion. Ademas, la cantidad de 

10 almacenamiento requerida es limitada al indicar, en forma codificada, la 
correccion focal que debe aplicarse al reloj de muestreo para obtener el 
siguiente foco. Mas aun, cada codigo de correccion focal puede ser compartido 
por multiples muestras consecutivas. Finalmente, la arquitectura de la tecnica 
aqui patentada presenta una elevada regularidad, facilitando su integracion en 

15 dispositivos VLSI y su realizacion modular. 

Se reclaman, como propias, las siguientes reivindicaciones: 

1. Un metodo de conformacion de haces en recepcion, denominado Tecnica 
de Correccion Focal Progresiva, para aplicaciones en las que una 
diversidad de N transductores, agrupados en un array, reciben energfa 

20 vibratoria desde focos F h i=0, 7, 2,..., situados a lo largo de una direccion 
programable, caracterizado porque muestrea la senal recibida en cada 
canal en el instante de llegada desde cada foco con un error absolute 
inferior a medio periodo de un reloj maestro de periodo T x , de forma que las 
muestras consecutivas e ki , i=0, 1, 2, ... obtenidas en cada canal k 

25 correspondan a los valores consecutivos de la senal originada en cada foco 
F\, i=0, 1, 2,..., con lo que cada muestra de la senal resultante n, i=0, 1, 2, ... 
se realiza mediante la suma de las N muestras e k i, k=1, 2, .., N : 

N 

r i = Y* e *i i = 0,l,2,... 

k=\ 

o de una combinacion lineal como: 



